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The  study  of  cathode  phenomena  in  gaseous  discharges  and 
electrical  arcs  has  progressed  over  the  last  hundred  years.  In 
the  last  twenty  years,  this  research  has  been  revived  as  the 
interest  in  high  power  devices  for  pulsed  power  applications, 
switching  devices,  and  nuclear  fusion  power  generation  has 
increased.  The  discovery  in  1958  of  a  new  form  of  electrical 
arc,  one  that  burns  between  a  metal  plate  and  a  dense  plasma, 
called  a  unipolar  arc  has  further  increased  interest  in  the 
cathode  phenomena.  For  both  the  bipolar  and  unipolar  arcs,  the 
current  transfer  occurs  in  localized  regions  between  a  cathode 
and  plasma  which  is  referred  to  as  the  cathode  spot.  It  has 
become  increasingly  apparent  that  the  phenomena  occurring  in 
the  cathode  spot  of  a  unipolar  arc  and  that  occurring  in  a 
bipolar  arc  are  similar.  Therefore,  much  of  the  theoretical  work 
done  on  the  cathode  spot  of  a  bipolar  arc  can  be  applied 
towards  the  understanding  and  formulation  of  a  consistent 
theory  of  the  unipolar  arc. 

This  study  will  begin  with  a  short  review  of  some  basic 
theories  of  plasma  physics,  space  charge  effects,  and  electron 
emission  mechanisms.  Two  useful  unipolar  arcing  models  are  then 
reviewed.  The  stationary  cathode  spot  will  then  be  studied  and 
a  set  of  equations  presented  that  can  be  used  to  determine  it's 


parameters.  It  will  be  found  that  the  system  of  equations  is 
not  closed,  however,  a  method  for  the  solution  of  the  system 
can  be  formulated  by  using  the  Steenbeck  Minimum  Principle. 
Finally,  results  of  calculations  based  on  these  equations  is 
presented. 


II. 


TH8QBY 


A.  PLASMA  PHYSICS 

L  Plasma  Characteristics 

A  plasma  is  formed  when  atoms  are  ionized  into  positive 
ions  and  negative  electrons  such  that  they  form  a  gas. 
However,  not  all  such  gases  are  considered  to  be  plasmas.  A 
more  complete  definition  of  a  plasma  has  been  given  by  Chen 
[Ref.  1]: 

"A  plasma  is  a  quasineutral  gas  of  charged  and  neutral 

particles  which  exhibits  collective  behavior.  " 

Collective  behavior''  means  that  the  behavior  of  a  small  region 
of  the  plasma  is  determined  not  only  by  conditions  in  the  plasma 
that  immediately  surrounds  it,  but  also  by  conditions  in  the 
plasma  far  away  from  it. 

A  plasma  can  be  characterized  by  two  quantities,  n,  the 
plasma  density  and  T,  the  plasma  temperature.  Since  the  plasma 
must  be  quasineutral  the  density  of  ions,  m,  and  the  density  of 
electrons,  ne,  must  be  approximately  equal  Therefore,  the 
plasma  density  is  defined  as 

nm  n,  *  nt.  (  eqn.  2.1  ) 

The  plasma  temperature  is  generally  given  in  units  of  energy 
which  corresponds  to  kT,  where  k  is  Boltzmann's  constant  and  T 
is  the  plasma  temperature  in  Kelvin.  For  a  gas  in  thermal 


11 
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equilibrium,  which  follows  a  Maxwellian  distribution,  the  average 
energy  of  the  gas  particles  is  related  to  the  plasma 
temperature  by 


E 


av 


C  eqn.  2.2  ) 


where  the  plasma  is  assumed  to  be  three  dimensional 
2.  Debye  ShtelH-ing 

When  a  plasma  is  subjected  to  an  external  electric  field 
the  electrons  and  ions  will  drift  in  opposite  directions.  As 
they  drift  apart,  they  produce  a  charge  separation  and  an 
internal  electric  field,  which  opposes  the  external  field.  The 
width  of  the  region,  over  which  the  charge  separation  must  occur 
to  balance  the  external  electric  field,  is  proportional  to  the 
thermal  energy  of  the  plasma  particles.  If  the  dimension  of  the 
bulk  plasma  is  greater  than  the  dimension  of  the  region  over 
which  this  charge  separation  occurs,  the  interior  of  the  plasma 
will  be  shielded  from  the  external  field. 

The  width  of  the  charge  separation  can  be  determined, 
approximately,  by  considering  a  plasma  where  the  ions  are  fixed 
in  space,  over  the  time  frame  of  interest,  and  the  electrons 
obey  a  Maxwellian  distribution  such  that  the  Boltzman  relation 
applies.  Using  Poisson's  Equation  in  one  dimension  gives 

d2v  e  r  s  (  e(in-  2-3  ) 

— - \n.-nA 
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where  V  is  the  electric  potential  and  e  is  the  electric  charge. 
Since  the  ions  are  considered  fixed,  the  density  of  ions  will  be 
constant  and  equal  to  the  plasma  density  no.  The  Boltzman 
relation  then  gives  the  electron  distribution  in  the  region  of 
interest  as 


(  eqn.  2.4  ) 


Substituting  equation  2.4  into  equation  2.3  gives 


d2V 

dx2 


-Wff  1- 

€o  \  \kTiJ 


(  eqn.  2.5  ) 


To  solve  equation  2.5  the  exponential  can  be  expanded  in  a 
Taylor  series,  and  neglecting  terms  of  second  order  and  higher 
gives 


d2V  fie2  (eqn.  2.6) 

dx2  €0kTa 

The  solution  of  equation  2.6,  which  is  a  homogenous  second  order 
linear  ordinary  differential  equation  is  given  by 

f  |xH  C  eqn.  2.7  ) 

^0exp[-UJ 

where 


(  eqn.  2.8  ) 


This  quantity  is  known  as  the  Debye  Length  and  is  a  measure  of 
the  thickness  of  the  shield  that  screens  the  bulk  plasma  from 


the  affect  of  external  fields.  This  screening  applies  as  long  as 
the  dimensions  of  the  bulk  plasma  is  much  greater  than  the 
Debye  Length. 

3.  Sheaths 

When  a  plasma  comes  in  contact  with  a  wall,  the  electrons 
and  ions  that  hit  the  wall  will  recombine.  But  since  the 
electrons  are  moving  at  a  higher  velocity  than  the  ions,  the 
electrons  will  be  lost  faster  than  the  ions  resulting  in  a  net 
positive  charge  at  the  plasma-wall  boundary.  The  wall  will  then 
be  at  a  lower  potential  than  the  bulk  plasma  and  an  electric 
field  will  exist.  Figure  2.1  illustrates  the  potential  variation  in 
a  plasma  which  is  in  contact  with  a  wall  Due  to  Debye  Shielding 
a  layer  of  charge  separation  will  exist,  called  the  sheath,  which 
will  isolate  the  bulk  plasma  from  this  potential  difference.  The 
effect  of  the  sheath  is  to  accelerate  the  ions  and  to 
decelerate  the  electrons  that  enter  the  region  until  the  flux  of 
ions  is  balanced  by  the  flux  of  electrons.  Therefore,  only  those 
electrons  in  the  high  energy  tail  of  the  velocity  distribution 
will  be  energetic  enough  to  cross  the  sheath  and  all  others  will 
be  repelled  back  into  the  plasma. 

The  variation  of  the  potential  in  the  plasma  sheath  will 
now  be  considered.  In  order  to  simplify  the  problem  the 
following  assumptions  will  be  made: 

*  The  ions  enter  the  sheath  region  with  a  drift  velocity 
uo , 

*  The  ions  have  Ti=0  so  that  all  ions  have  a  velocity  uo 
at  the  plasma- sheath  boundary, 


*  The  sheath  region  is  collisionless  and  in  steady  state, 

*  The  potential  decreases  mo  no  tonic  ally  with  x, 

*  The  electrons  follow  the  Boltzman  relation  (equation 
2.4). 


Figure  2.1  Plasma-Wall  Interaction 


Applying  conservation  of  energy  to  the  ions  in  the  sheath 
region  gives 

^A/u2-  il\-9V{x).  (eqn.2.9) 

t 

Solving  equation  2.9  for  the  ion  velocity  gives  1 


u 


(  eqn.  2.10  ) 


The  ion  density  within  the  sheath  is  determined  from  the  ion 
equation  of  continuity  which  is  given  by 


1 


Using  equation  2.10  for  the  ion  velocity,  equation  2.11  can  be 
solved  for  the  ion  density  which  gives 

-i  (  eqn .  2.12  ) 

2eV\ 


nt(x) 


n0(l- 


Uo  M  ) 

Substituting  equation  2.4  and  2.12  into  equation  2.3  gives 


d2V  eno 

€0 


dx‘ 


exp 


qV 


kT 


i)-(-g) 


(  eqn.  2.13  ) 


Equation  2.13  is  the  non-linear  planar  sheath  equation.  By 
multiplying  each  side  by  the  first  derivative  of  the  potential 
the  equation  can  be  integrated  once  to  give 


This  inequality  can  be  greatly  simplified  by  expanding  the  left 
hand  terms  in  Taylor  series  expansions  and  neglecting  terms  of 
third  order  and  higher.  The  resulting  inequality  is  known  as  the 
Bohm  sheath  criterion  and  is  given  by 

(  eqn.  2.16  ) 


„  fkT. 


This  requires  that  the  ions  must  enter  the  sheath  region  with  a 
velocity  greater  than  the  ion  acoustic  velocity  of  the  plasma. 
In  order  for  this  to  occur,  there  must  be  a  finite  electric  field 
within  the  plasma.  Therefore,  the  assumption  that  the  potential 
and  the  electric  field  are  zero  at  the  plasma-sheath  boundary  is 
only  an  approximation. 


B.  SPECS  CHARGE  LIMITED  CURRENT  FLOW 

1  Child— Lamnwiir  Law 

The  current  that  can  flow  between  the  cathode  and  anode 


of  a  vacuum  diode  is  limited  by  the  space  charge  of  the 
electrons  that  exist  between  them.  The  electrons  distort  the 
external  field  and  thereby  reduce  and  can  even  reverse  the  field 
at  the  cathode  surface. 

As  an  example,  consider  a  diode  consisting  of  infinite  flat 
plates  separated  by  a  distance  d.  Assume  that  there  is  an 
unlimited  supply  of  electrons  available  from  the  cathode,  and 
that  the  initial  velocity  of  the  electrons  after  emission  is 
negligible  compared  with  the  velocity  that  they  will  gain  while 


crossing  the  electrode  gap.  The  kinetic  energy  of  the  electrons 
crossing  the  gap  and  the  current  density  carried  by  the 
electrons  is  given  by 


^rnu2  -  aV 


j  -  env. 


(  eqn.  2.17  ) 


(  eqn.  2.18  ) 


Using  Poisson's  equation,  in  one  dimension,  the  current- voltage 


relationship  can  be  derived  as  follows: 

d2 V (x)  _  an (x) 
dx2  e0 


(  eqn.  2.19  ) 


The  electron  velocity  is  determined  from  equation  2.17  to  be 


v(x) 


•yfe  v'\ 

V  m 


(  eqn.  2.20  ) 


The  electron  density  is  determined  using  equation  2.18  and  2.20, 


and  after  substituting  into  equation  2.19  gives 


v.l  rEv-1* 

:2  e0V2fl  * 


(  eqn.  2.21  ) 


In  order  to  solve  equation  2.21  both  sides  are  multiplied  by  the 
first  derivative  of  the  potential  and  reduced  to  the  following 


form 


(dvV  i  [2mv-idV 

;i3ij  r,y  ~v  di- 


(  eqn.  2.22  ) 
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Equation  2.22  is  then  integrated  over  the  gap  using  the 
boundary  conditions  that  the  electric  field  and  the  potential  are 
zero  at  the  cathode.  The  result  of  this  integration  gives 


(  eqn .  2.23  ) 


Integrating  once  again  over  the  gap,  applying  the  boundary 
conditions,  and  rearranging  gives  the  desired  current- voltage 
relationship  of 


(  eqn.  2.24  ) 


Equation  2.24  is  known  as  the  Child- Langmuir  Law  for 
space  charge  limited  current  flow  [Ref.  2].  Although  it  has  been 
derived  here  using  simple  assumptions,  the  proportionality  of  the 
current  density  to  the  3/2  power  of  the  potential  difference 
remains  for  more  difficult  geometries  under  non- relativistic 
conditions.  It  should  also  be  pointed  out  that  the  Child- Lang¬ 
muir  Law  applies  to  both  positive  and  negative  charge  can 
Figure  2.2  is  a  plot  of  electron  current  density  ver* 
potential  for  several  gap  widths. 

2.  Space.  Charge  Llalted  Current  in  a.  Sheath 

In  deriving  the  planar  sheath  equation  it  was  assumed 
that  the  electron  density  followed  the  Boltzroan  relation 
equation  2.4.  Therefore,  in  a  region  close  to  the  wall  the 
electron  density  will  be  much  less  than  the  ion  density 


3.  jUriK^nn  Ba.mt.Jrm 

In  an  arc,  the  current  is  carried  primarily  by  the 

electrons  emitted  by  the  cathode  and  by  a  smaller  amount  of 
positive  ions  flowing  from  the  plasma  to  the  cathode.  If  the 
space  charge,  which  causes  the  cathode  drop  in  the  plasma 
sheath,  is  high  enough  it  may  exert  a  strong  electric  field  at 

the  surface  of  the  cathode.  This  can  then  result  in  a  large 

field  emission  current  in  addition  to  any  thermionic  current. 

The  electric  field  at  the  cathode  surface  can  be 

determined  using  Poisson's  Equation.  The  following  derivation 
measures  distance  and  potential  difference  from  inside  the 
plasma  where  the  cathode  drop  begins.  The  electron  and  ion 
current  densities  are  given  by 

j %  *  Qfl'V  C  eqn.  2.29  ) 

and 

ji  *  ontu .  (  e<ln-  2-30  ) 

Substituting  equation  2.29  and  2.30  into  equation  2.3  gives 

d2V _ if  (  eqn<  2,31  > 

dx2  €q \U  v) 

where  V  is  the  potential  difference,  ji,  the  ion  current  density, 
je,  the  electron  current  density,  u,  the  ion  velocity,  and  v,  the 
electron  velocity. 

The  velocity  of  the  ions  and  electrons  are  determined 
using  the  conservation  of  energy.  The  initial  energies  are 


neglected,  based  on  the  assumption  that  the  energy  gained  by 
the  ions  and  electrons  in  crossing  the  sheath  are  much  greater 
than  the  initial  energies.  For  the  ions  the  energy  balance  is 


A/  12  2  =  eV 

2 

and  for  the  electrons  the  energy  balance  gives 

^mv2  =  e(vc- V) 


(  eqn.  2.32  ) 


(  eqn.  2.33  ) 


where  M  is  the  mass  of  the  positive  ion,  m,  the  mass  of  the 
electron  and  Vc  is  the  potential  of  the  cathode.  Substituting 
equation  2.32  and  2.33  into  equation  2.31  gives 


d2V 

dx2 


.  (  M  V  f  m 
~e0  J\2eV)  ~ J\2e{Vc - 


(  eqn.  2.34  ) 


Multiplying  both  sides  by  the  first  derivative  of  the  potential 
difference,  equation  2.30  can  be  integrated  once  and  gives 


-2  4)  .(UV\*  .  f  m(K«-K)V  .(mV, 
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(  eqn.  2.35  ) 

where  F  is  the  electric  field.  Evaluating  for  V  -  Vc  gives  the 


electric  field  at  the  cathode  surface  as 

.  4(UV,\i  (  i.fm'SX 


(  eqn .  2.36  ) 


2»' 

y 

ft*: 


Equation  2.36  is  known  as  Macheown's  Equation  [Ref.  3]  and 
can  be  used  to  determine  the  electric  field  existing  at  the 
cathode  provided  the  values  of  the  cathode  drop,  the  ion 
current  density,  and  the  electron  current  density  are  known. 
Figure  2.3  is  a  plot  of  electric  field  at  the  cathode  surface 
versus  ion  current  density  for  several  ratios  of  electron  to 
ion  current  densities. 


C.  ELECTRON  EMISSION  FROM  METALS 

The  theory  of  electron  emission  from  metals  is  very  well 
understood.  The  basic  processes  have  been  explained  on  the 
basis  of  a  free  electron  gas  in  the  metal,  that  is  described  by 
Fermi- Dirac  statistics,  and  a  potential  barrier  at  the  metal 
surface.  Figure  2.4  shows  a  typical  potential  energy  diagram 
for  an  electron  as  a  function  of  distance  from  the  metal- vacuum 
interface.  For  thermionic  emission,  that  is  emission  at  high 
temperature  and  low  electric  field,  the  electrons  are  emitted 
predominately  by  passing  over  the  potential  barrier  at  the  metal 
surface.  For  field  emission,  that  is  emission  at  low  temperature 
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Figure  2.3  MacKeown's  Equation  for  Electric  Field  at  the 

Cathode  Surface 
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Figure  2.5  Regions  of  Applicability  for  Thermionic,  Field,  and 
Intermediate  Electron  Emission  Expressions 

For  equation  2.41  to  be  applicable  the  following  conditions  must 

be  satisfied: 


(  eqn.  2.43  ) 
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(  eqn.  2.44  ) 


d  I  (1-d) 

Equation  2.41  reduces  to  the  Richardson  equation  in  the  limit  as 
the  electric  field  is  reduced  toward  zero.  Appendix  D  is  a 
listing  of  a  Fortran  subroutine  that  calculates  the 
Thermionic -Fie  Id  emission  current  density  based  on  equation  2.4L 

2.  Field- Thermion-in  Raisa  ion 

The  current  density  that  can  be  expected  from  field 


1 


Equations  2.46,  2.48,  and  2.49  are  the  Fowler- Nordheim  Elliptic 
functions.  The  functions  K[k]  and  E[k]  are  complete  elliptic 
integrals  of  the  first  and  second  kind  respectively  and  are 
defined  as 


m 

KW-ftl  -k2s\n2e)'idO, 


(  eqn.  2.51  ) 


m 

E[k]~  J*(  1  -  fc2sin 2Q)* dLQ. 


(  eqn.  2.52  ) 


For  equation  2.45  to  be  applicable  the  following  conditions  must 
be  satisfied: 


J  F3/*  kT 

<t-  F  > - + - , 

v  it  1  -ckT 


(  eqn.  2.S3  ) 


\-ckT>f2fkT, 
fm  2 F{t2-F)V^y)' 


(  eqn.  2.54  ) 


(  eqn.  2.S5  ) 


Equation  2.45  reduces  to  the  Fowler- Nordheim  equation  in  the 
limit  as  the  temperature  approaches  zero.  Appendix  E  is  a 
listing  of  a  Fortran  subroutine  that  calculates  the 
Fie  Id- Thermionic  current  density  based  on  equation  2.45.  This 
subroutine  makes  calls  to  the  subroutines  listed  in  Appendix  F, 
G,  and  H,  which  are  used  to  calculate  the  Complete  Elliptic 
Integral  of  the  First  Kind,  Complete  Elliptic  Integral  of  the 
Second  Kind,  and  the  Fowler-Nordheim  Elliptic  functions 
respectively. 


Appendix  I  is  a  listing  of  a  Fortran  subroutine  that 
calculates  the  intermediate  Fie  Id- Thermionic  emission  current 


density  based  on  equation  2.56. 

4.  BxplagiVB  Electron  Bilaalon 

The  emission  mechanisms  discussed  so  far  have  wide 
applicability  to  many  physical  devices.  However,  with  the  advent 
of  high  current  cathodes  the  thermionic  and  field  emission 
theories  were  unable  to  explain  the  high  current  densities 
observed.  One  explanation  for  this  high  current  density  is  that 
the  surface  of  the  cathode  is  not  smooth,  but  actually  has  a 
surface  covered  with  micropoints,  or  whiskers.  The  significance 
of  the  whisker  is  that  the  electric  field,  at  the  surface  of  the 
whisker  is  enhanced  by  a  factor  of  10  -  1000.  Figure  2.6 
illustrates  the  electric  field  enhancement  that  can  occur  at  the 
tip  of  a  whisker.  This  results  in  a  higher  field  emission 
current  density  from  the  whisker,  and  ultimately  leads  to 
explosive  destruction  of  the  emission  site.  This  destruction 
produces  a  local  burst  of  plasma  that  expands  rapidly  away  from 
the  cathode  surface  and  has  therefore  come  to  be  known  as  a 
cathode  flare  or  plasma  jet.  The  electrons,  that  reach  the 
anode,  are  then  emitted  from  the  surface  of  the  plasma  jet. 
This  whole  process  is  commonly  referred  to  as  explosive  electron 
emission  [Refs.  5  and  6]. 
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Figure  2.6  Electric  Field  Enhancement  at  a  Whisker 


Fursei  and  Zhukov  [Fefs.  7,  8,  and  9]  have  experimentally 
characterized  the  explosive  emission  process  and  have  confirmed 
the  existence  of  raicropoints  and  that  the  current  density 
follows  the  Fowler- Nordheim  field  emission  theory  with  the 
inclusion  of  a  field  enhancement  factor. 


In  1958,  Robson  and  Thoneman  observed  an  electrical  arc 
burning  between  a  dense  hot  plasma  and  a  metal  wall  [Ref.  10]. 
What  made  this  unusual  was  that  there  was  no  anode  present  to 
collect  the  electrons  emitted  from  the  conducting  wall 
Apparently,  the  wall  was  acting  as  both  a  cathode  and  anode,  and 
therefore,  they  named  what  they  observed  a  ‘unipolar  arc." 

Robson  and  Thoneman  proposed  the  following  model  to  explain 
the  occurrence  of  unipolar  arcing.  When  a  hot  dense  plasma 
comes  in  contact  with  a  conducting  wall  a  sheath  will  form.  If 
the  plasma  temperature  is  high,  the  sheath  or  floating  potential 
can  exceed  that  necessary  to  initiate  an  electric  arc.  A 
cathode  spot,  which  is  a  local  center  of  emission,  could  then 
form  on  the  metal  surface.  The  large  electron  current  flow  into 
the  plasma  from  the  arc  would  result  in  lowering  the  plasma 
potential  to  the  cathode  fall  potential  of  the  arc.  This  allows 
more  electrons  from  the  high  energy  tail  of  the  velocity 
distribution  to  cross  the  sheath  and  recombine  at  the  wall  and 
thus  closing  the  current  loop.  Figure  3.1  illustrates  the 
situation  in  the  sheath  prior  to  initiation  of  an  arc,  and  Figure 
3.2  shows  the  return  current  flow  after  arc  initiation. 
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where  I  is  the  current,  A  is  the  area  of  the  metal  plate,  and  Vc 
is  the  cathode  drop  of  the  arc.  But  since  the  net  current  flow 
between  the  plasma  and  the  wall  must  be  zero,  equation  3.1  also 
gives  the  current  flow  inside  the  arc. 

The  ion  current  density  in  the  sheath  can  be  approximated 


Ji~en  o 


rkTa 


(  eqn .  3.2  ) 


where  the  Bohm  sheath  criterion  is  used  to  approximate  the  ion 
velocity.  The  electron  current  density  in  the  sheath  can  be 
determined  by  noting  that  the  net  flux  of  electrons,  *,  from  the 
plasma  into  the  sheath  is  given  by 


(  eqn.  3.3  ) 


Using  the  Boltzman  relation  for  the  electron  density  and  the 
expression  for  the  average  speed  of  particles  following  a 
Maxwellian  distribution  gives 


y.-5<m0exp^  — 


QkT , 


(  eqn.  3.4  ) 


The  net  current  density  across  the  sheath  must  be  zero. 
Therefore,  setting  equation  3.2  equal  to  equation  3.4,  and 


solving  for  the  sheath  potential  gives 

jq. 

f  2e  \2nmJ 


(  eqn.  3.5  ) 


Robson  and  Thoneman  point  out  that  there  is  a  minimum 
current  necessary  to  maintain  an  arc  burning.  Therefore, 
according  to  equation  3.1,  if  the  plate  area  is  small  enough  a 
unipolar  arc  can  not  occur.  For  example,  using  copper  with  a 
plasma  temperature  of  4  eV,  a  cathode  fall  potential  of  15  V,  a 
plasma  density  of  1013  cm- 3 ,  and  an  arc  current  of  100  amperes 
would  require  a  plate  area  of  115  cm2.  For  a  circular  area,  the 
radius  required  is  6.05  cm. 

B.  Schwirzke-Taylor  Model 

Schwirzke  and  Taylor  have  expanded  on  the  basic  concepts  of 
the  Robson- Thoneman  model  [Ref.  11].  They  propose  that  for  an 
arc  to  occur,  the  ion  density  above  the  cathode  spot  must 
increase  to  allow  for  the  arc  current  to  flow.  This  is  contrary 
to  the  Robs  on- Thoneman  model  which  assumes  a  constant  plasma 
density.  Once  a  cathode  spot  forms  the  desorption  of  neutral 
atoms  and  evaporation  of  the  metal  atoms  from  the  spot  is 
greatly  increased.  If  only  a  small  fraction  of  these  neutral 
atoms  become  ionized  in  the  sheath  the  plasma  density  will 
increase  above  the  cathode  spot.  Thus,  according  to  equation 
2.8  the  Debye  length  will  decrease  resulting  in  a  decrease  of 
the  sheath  width  above  the  cathode  spot.  The  electric  field 
above  the  spot  increases,  where  the  field  can  be  approximated 


Schwirzke  [Refs.  12  and  13]  elaborates  further  on  the 
electric  fields  set  up  inside  the  high  density  plasma  which  forms 
above  the  spot.  The  fluid  equation  of  motion  of  the  electrons 


in  the  plasma  is  given  by 


<?vt  /-  -  —x 


(  eqn .  3.7  ) 


The  time  derivative  can  be  neglected  over  the  time  span  of  an 
arc,  and  further  assuming  no  magnetic  field,  equation  3.7  can  be 
solved  for  the  electric  field,  which  gives 


p  =  Z_  J_v/> 

a  ona 


(  eqn.  3.8  ) 


The  first  term  on  the  right  hand  side  of  equation  3.8  can  be 
neglected  in  a  weakly  ionized  plasma.  This  implies  that  the  mean 
free  path  length  for  collision  between  the  electrons  and  the 
neutral  atoms  is  much  greater  than  the  dimensions  of  the  region 
of  increased  density  above  the  cathode  spot.  Thus,  the 
increased  plasma  pressure  above  the  spot  creates  a  radial 


electric  field  given  by 

1  kT , 

F - V/># - -Vn8. 

en,  ena 


(  eqn.  3.9  ) 


The  affect  of  this  radial  electric  field  is  to  lower  the  plasma 
potential  in  a  ring  around  the  cathode  spot.  The  change  in  the 
potential  can  be  determined  by  integrating  equation  3.9  from  the 
center  of  the  arc,  where  the  plasma  density  is  highest,  radially 
outward  until  the  plasma  density  equals  the  background  density. 
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The  change  in  potential  is  then  given  by 

e  \noQJ 


(  eqn.  3.10  ) 


In  this  ring  of  reduced  potential,  electrons  in  the  tail  of  the 
Maxwellian  distribution  can  cross  the  sheath  more  easily  and 
thus  a  higher  return  current  flow  from  the  plasma  to  the  wall 
results.  This  current  flow  closes  the  current  loop  of  the 
unipolar  arc.  The  concepts  of  the  Schwirzke- Taylor  model  are 
illustrated  in  Figure  3.3. 
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Figure  3.3  Schwirzke- Taylor  Unipolar  Arc  Model 
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If  equation  3.10  is  set  equal  to  the  floating  sheath 
potential,  equation  3.5,  the  ratio  of  the  electron  density  in  the 
plasma  above  the  spot  to  the  background  plasma  density  which 
reduces  the  sheath  potential  to  zero  can  be  determined  as 


(  eqn.  3.1  1  ) 


n90  v  2jrm 


Equation  3.3  and  equation  3.11  can  be  used  to  calculate  the 
maximum  area  necessary  to  support  a  given  unipolar  arc  return 
current.  By  setting  the  current  equal  to  the  current  derived 
for  the  Robson-Thoneman  model  equation  3.1,  the  ratio  of  the 
return  current  areas  required  by  each  model  can  be  compared. 


The  result  of  the  comparison  gives 


U-r  /  M  /  (  elO  [2nm\  1  /  hi 

U-r  *  V  2irm\CXPV  kTJ  \  hi  f  >  \  2 nm 


(  eqn.  3.12  ) 


where  Ar-t  and  As-T  are  the  areas  required  by  the 
Robson-Thoneman  model  and  Schwirzke-Taylor  model  respectively. 
As  an  example,  for  copper  with  a  cathode  fall  potential  of  15 
volts  and  a  plasma  temperature  of  4  eV  the  area  ratio  is  8.43  x 
103.  Therefore,  for  a  100  ampere  arc  the  area  required  by  the 
Schwirzke-Taylor  model  is  only  L37  x  10*2  Cm2,  or  for  a  circular 
area  the  radius  required  is  only  6.59  x  10*2  cm. 
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IV.  CATHODE  SPOTS 


A.  TYPES  OF  CATHODE  SPOTS 

A  cathode  spot  is  a  small,  luminous  region  on  the  cathode, 
where  a  localized  transfer  of  current  occurs  between  the 
cathode  and  a  plasma  in  an  arc  discharge.  The  type  of  arc 
discharge  can  be  either  bipolar  or  unipolar,  however  in  all  cases 
a  cathode  spot  is  present.  In  addition,  the  spots  are 
associated  with  the  erosion  of  the  cathode  and  therefore,  are 
also  localized  centers  for  the  transfer  of  material  from  the 
cathode  into  the  plasma.  This  has  been  evidenced  by  the 
formation  of  craters  on  the  surface  of  the  cathode. 

It  has  been  possible  over  many  observations  to  distinguish 
between  two  different  types  of  cathode  spots  [Ref.  14].  The 
first  type,  or  type  I  spots,  are  characterized  as  rapidly  moving 
spots  with  speeds  of  approximately  10«  cm/s  across  the  surface 
of  the  cathode  and  with  relatively  low  erosion  from  the  surface. 
The  type  II  spots,  on  the  other  hand,  have  speeds  of 
approximately  102  cra/s  and  have  much  greater  erosion  from  the 
surface.  The  type  II  spot  also  has  a  tendency  to  group  together 
and  form  clusters.  The  mechanism  of  erosion  for  the  type  I 
spot  is  believed  to  be  due  to  non-thermal  means,  such  as 
explosions  of  whiskers  or  micro  inhomogenie  ties  on  the  cathode 
surface.  The  erosion  mechanism  for  the  type  EE  spot  is  believed 


to  be  due  to  thermal  evaporation. 


B.  THKOBT  OF  STATIONABY  CATHODE  SPOTS 

The  theoretical  modeling  of  cathode  spots  has  progressed 
over  the  last  100  years,  but  up  until  the  present  time  there 
has  been  no  satisfactory  understanding  of  the  processes 
occurring  within  the  spot.  The  difficulty  occurs  due  to  the 
diverse  phenomena,  which  in  themselves  are  extremely  complicated 
to  describe  mathematically,  that  take  place  within  the  small 
cathode- plasma  region  known  as  a  cathode  spot.  A  detailed 
analysis  would  require  an  extensive  use  of  the  kinetic  theories 
of  Plasma  Physics,  Solid  State  Physics,  Physics  of  High 
Temperature  Phenomena,  and  Gas  Dynamics.  In  order  to  avoid 
some  of  the  complications  the  majority  of  models  developed  for 
the  cathode  spot  are  stationary,  and  therefore  are  unable  to 
explain  the  initiation,  motion,  and  disappearance  of  the  cathode 
spot  [Refs.  12  and  15]. 

The  energy  balance  method  developed  by  Lee  and  Greenwood 
[Ref.  16]  has  been  the  basis  of  the  majority  of  stationary 
models  including  those  proposed  by  Kulyapin  [Ref.  17],  Kubono 
[Ref.  18],  Beilis  [Ref.  19],  and  Moizhes  and  Nemchinskii  [Refs.  20, 
21,  22,  and  23].  The  latter  work  by  Moizhes  and  Nemchinskii  are 
the  basis  for  the  model  developed  in  this  paper.  The  cathode 
spot  is  divided  into  three  regions  as  illustrated  in  Figure  4.L 
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Figure  4.1  Model  of  the  Cathode -Plasma  Region 


Region  I,  referred  to  as  the  Langmuir  region,  is  a  thin 
space  charge  sheath  over  which  the  cathode  potential  drop,  Vk, 
occurs.  It  is  in  this  region  that  the  electrons  emitted  by  the 
cathode  are  accelerated  into  the  plasma  and  where  the  ions 
produced  in  the  generation  region  are  accelerated  towards  the 
cathode.  Region  II,  referred  to  as  the  generation  region  or 

Knudsen  Region,  is  where  the  neutral  atoms  evaporated  from  the 
cathode  are  ionized.  Region  HI  is  an  expanding  channel  or 

transition  region  over  which  the  highly  dense  current  channel  in 
regions  I  and  II  expands  out  into  the  bulk  plasma. 
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In  developing  a  theoretical  model  of  a  cathode  spot,  one 
of  the  first  decisions  that  must  be  made  is  the  electron 
emission  mechanism.  In  general  most  models  use  a  combination 
of  thermal  and  field  equations  that  are  expressed  as 

j.-f{T,.F)  (eqn.4.1) 

where  je  is  the  electron  current  density,  Tk  is  the  cathode 
temperature,  and  F  is  the  electric  field  at  the  cathode  surface. 
In  actual  use,  expressions  such  as  equation  2.41,  2.45,  or  2.56 
are  used.  The  particular  equation  used  depends  on  the 
applicability  under  the  given  conditions.  The  electric  field  at 
the  cathode  surface  is  determined  using  MacKeown's  Equation, 
equation  2.36. 

An  energy  balance  must  exist  in  the  plasma  region  above 
the  cathode  spot  and  can  be  described  by  equation  4.2. 


(  eqn.  4.2  ) 


Energy  is  brought  into  the  region  by  the  electrons  emitted  from 
the  cathode  and  consists  of  the  thermal  energy  and  the  energy 
gained  as  the  electrons  were  accelerated  through  the  cathode 
potential  drop.  Energy  is  removed  from  the  region  by  ions 
flowing  back  to  the  cathode.  This  energy  loss  consists  of  the 
ionization  energy  and  the  thermal  energy  of  the  ions.  A  further 
loss  occurs  from  both  ions  and  electrons  that  leave  the  region 
and  pass  into  the  bulk  plasma  and  is  described  by  the  Peltier 


heat  flux  jfi kT,.  In  a  plasma,  in  which  Coulomb  collisions  are  the 
predominant  scattering  process,  the  energy  transfer  coefficient, 
f,  is  approximately  equal  to  3.21  It  should  be  noted  that 
electrons  in  the  plasma,  also,  can  leave  the  region  towards  the 
cathode,  but  only  if  energetic  enough  to  over  come  the  cathode 
potential  drop,  Vk.  This  reverse  electron  current  can  be 
approximated  by 

QTjaP  f  (  eqn.  4.3  ) 

y'“~V2,r*rpeXPr^J- 

For  low  plasma  temperatures  and  high  cathode  potential  drops 
the  reverse  electron  current  density  is  much  less  than  the  ion 
current  density  and  can  therefore  be  neglected  in  this  model 

An  energy  balance  at  the  cathode  surface  must  also  exist 
and  can  be  described  by  equation  4.4. 

(  eqn.  4.4  ) 

The  predominant  loss  mechanisms  considered  are  heat  conduction, 
emissive  cooling,  and  evaporation.  The  primary  energy  source 
considered  is  ion  bombardment  of  the  cathode  surface.  The  ions 
bring  to  the  surface  thermal  energy,  ionization  energy,  and 
energy  gained  while  being  accelerated  by  the  cathode  potential 
drop.  Energy  is  released,  approximately  equal  to  the  work 
function,  when  ions  and  electrons  recombine  at  the  surface.  The 
electrons  emitted  by  the  surface  remove  thermal  energy  and  an 


4KRTk  ■  JtR2 
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energy  approximately  equal  to  the  work  function.  In  addition, 
neutral  atoms,  evaporating  from  the  surface,  remove  energy 
equivalent  to  the  heat  of  vaporization.  Heat  conduction  from 
the  cathode  spot  is  based  on  a  circular  spot  of  radius  R  at  a 
constant  temperature  Tk  and  where  k  is  the  thermal  conductivity 
of  the  cathode.  The  work  function  to  be  used  in  equation  4.4  is 
corrected  for  the  Schottky  effect  and  is  given  by 


*-*o- 


(  eqn.  4.5  ) 


The  total  current  flow  in  the  arc  is  given  by 


/  -  jrrrt2(  jt  +  /<).  (  «qn-  4-6  ) 

A  fraction  of  the  neutral  atoms  evaporated  from  the  cathode 
will  return  to  the  cathode.  The  return  coefficient  can  be 
determined  from  equation  4.7  which  is  based  on  a  Monte  Carlo 
calculation  performed  by  Nemchinskii  [Ref.  241. 


The  ion  current  density  is  determined  from 


(  eqn.  4.6  ) 


1 

Ji  *  -6fltU 


(  eqn.  4.8  ) 


where 


(  eqn.  4.9  ) 


Substituting  equations  4.9  and  4.10  into  equation  4.8  gives  the 
ion  current  density  as 


qP 


(  eqn .  4. 1  1  ) 


‘  UnhdkTt 

where  P  is  the  vapor  pressure  of  the  cathode  material  and  «  is 
the  degree  of  ionization.  The  degree  of  ionization  is  determined 
from  Saha's  Equation  which  is  expressed  as 


a  - 


-i 


kT 


(  eqn.  4.12  ) 


The  erosion  rate  of  neutral  atoms  is  given  by 

g~f2 TukT}'  n)' 


(  eqn.  4.13  ) 


The  potential  drop  inside  the  quasineutral  plasma  of  the 
transition  region  is  given  by 


C  eqn.  4.14  ) 


where  z  is  the  average  charge  of  the  ions  in  the  region,  j  is 
the  total  current  density  given  by 

y-yt+y'. 

and  a  is  the  electric  conductivity  given  by 


(  eqn.  4.15  ) 


1  .53  x  10**! 
• - rill  7' 


(  eqn.  4.16  ) 


where  the  Coulomb  Logarithm  is  approximated  as 


In  A  -  23 


The  total  potential  drop  across  the  arc  is  given  by 

Vvk-Vk  +  U. 

2.  Method  of  Solution 


(  eqn.  4.17  ) 


(  eqn.  4.18  ) 


The  system  of  equations  presented  in  the  previous  section 
is  not  closed.  Therefore,  more  information  is  required  in  order 
to  solve  the  system.  This  will  be  accomplished  here  using  the 
Steenbeck  Minimum  Principle  [Refs.  25  and  26].  Simply  stated  this 
requires  that  the  total  potential  drop  across  the  arc  must  be  a 
minimum  for  a  stationary  arc  to  exist.  Given  an  arc  current  and 
cathode  temperature  the  system  of  equations  is  solved  for  the 
cathode  potential  drop,  plasma  temperature,  ion  current  density, 
electron  current  density,  erosion  rate,  potential  drop  in  the 
transition  region,  and  total  arc  potential  drop.  The  cathode 
temperature  that  defines  the  stationary  cathode  spot,  for  the 
given  arc  current,  is  that  cathode  temperature  that  minimizes 
the  total  potential  drop.  If  a  minimum  total  arc  potential  can 
not  be  found  for  a  given  arc  current,  then  a  stable  arc  can  not 
exist  for  this  current. 


|t 


The  first  step  in  the  solution  is  to  determine  the  plasma 
temperature,  Tp,  and  cathode  drop,  Vk.  To  accomplish  this  a 
guess  is  made  for  the  plasma  temperature.  The  partial  pressure 
of  the  cathode  material  is  determined  using  [Ref.  27] 

log  P  ■  rir-+  Slog  T,  ♦  CT t  +  D  (eqn.4.19) 

where  for  copper  A=  -1752x10*,  B=  -1.21,  C=  0,  and  D=  13.21.  The 
degree  of  ionization  is  determined  using  equation  4.12  and  the 
return  coefficient  is  determined  using  equation  4.6.  The  ion 
current  density  can  now  be  determined  using  equation  4.11 
Experimentally,  it  has  been  observed  that  the  cathode  drop  is 
approximately  15  volts,  so  this  will  be  used  as  a  first 
approximation  in  order  to  determine  the  electric  field  at  the 
cathode  surface  using  equation  2.36.  The  Schottky  corrected 
work  function  is  determined  using  equation  4.5  and  the  electron 


current  density  is  determined  using  equation  4.1  The  spot 


v“-7Ti 

(  eqn.  4.21  ) 

The  potential  is  then  determined  using  both  equation  4.20  and 
4. 21.  The  results  are  compared  and  the  plasma  temperature  is 
then  adjusted.  This  procedure  is  then  repeated  until  the 
potentials  determined  from  both  equations  are  in  agreement.  The 
conductivity  is  determined  using  equation  4.16  and  4.17.  The 
potential  drop  in  the  transition  region  is  then  determined  using 
equation  4.14.  Finally,  the  total  arc  potential  is  determined 
from  equation  4.18.  Appendix  J  is  a  Fortran  program,  which  uses 
the  method  described  above,  to  calculate  the  cathode  spot 
parameters. 
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The  Fortran  program  listed  in  Appendix  J  was  used  to 
calculate  the  parameters  of  the  cathode  spot  for  arc  currents 
of  100,  150,  175,  200,  300  and  400  Amperes.  The  program 
required  knowledge  of  several  material  parameters  which  are 
specified  in  Table  L 


TABLE  1 
CATHODE  SPOT  MODEL 

CONSTANTS 

Constant 

Value 

Copper 


63.54 


300.3 


7.726 


Material 


Ion  Mass  (gr/mole) 


Heat  of  Vaporization  (kilojoules) 


Ionization  Energy  (eV) 


Work  Function  (eV) 


Thermal  Conductivity  (W/m-K) 


Figure  5.1  illustrates  the  variation  of  the  total  arc 
potential  with  cathode  temperature  for  different  arc  currents. 
The  cathode  temperature  for  a  stable  cathode  spot,  and  thus 
the  parameters  of  the  cathode  spot,  is  determined  from  the 
minimum  of  the  total  arc  potential  It  can  be  seen  from  Figure 
5.1  that  a  stable  arc  does  not  occur  for  arc  currents  below  200 
Amperes. 


The  results  of  the  calculations  are  listed  in  Table  2  through 
7  for  arc  currents  of  100,  150,  175,  200,  300,  and  400  Amperes 
respectively.  In  addition  Figures  5.2  through  5.7  illustrate  the 
variation  of  the  cathode  potential,  transition  region  potential, 
and  total  arc  potential  for  each  arc  current. 

The  cathode  spot  parameters  for  stable  arc  currents  of  200, 
300,  and  400  Amperes  are  summarized  in  Table  8.  The  results 
are  in  general,  order  of  magnitude,  agreement  with  experimental 
data  [Ref.  23],  however  stability  has  been  indicated  at  arc 
currents  between  100  and  200  Amperes.  This  discrepancy  occurs 
since  the  material  constants  listed  in  Table  1  are  not  really 
constants  at  these  temperatures.  A  more  detailed  comparison  is 
not  warranted,  since  the  published  experimental  data  for  the 
cathode  spot  parameters  vary  widely.  Confirmation  of  the  models 
of  the  cathode  phenomena  is  greatly  limited  by  this  lack  of 
consistent  experimental  results. 


TABLE  2 

CATHODE  SPOT  PARAMETERS  FOR  A  100  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode  Plasma  Cathode  Transition  Total  Arc 

Temperature  Temperature  Potential  Potential  Potential 

(Kelvin)  (eV>  (Volts)  (Volts)  (Volts) 


3500 

3550 

3600 

3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Cathode 
T  emperature 
(Kelvin) 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


22.89 

21.54 

20.34 
19.27 
18.32 
17.49 
16.75 
16.10 

15.54 
15.06 
14.66 

14.34 
14.09 
13.91 
13.82 


10.18 

10.65 

11.10 

11.53 
11.92 
12.26 

12.53 
12.71 
12.78 
12.70 
12.44 
11.94 


Partial 

Pressure 

(atm.) 

Degree  of 
Ionization 

Return 

Coefficient 

10.9 

0.1625 

0.4726 

12.6 

0.1712 

0.4751 

14.5 

0.1803 

0.4777 

16.6 

0.1897 

0.4804 

19.0 

0.1996 

0.4831 

21.6 

0.2101 

0.4859 

24.5 

0.2213 

0.4888 

27.6 

0.2332 

0.4918 

31.1 

0.2461 

0.4950 

34.9 

0.2602 

0.4983 

39.1 

0.2758 

0.5018 

43.6 

0.2931 

0.5057 

48.5 

0.3125 

0.5098 

53.8 

0.3348 

0.5144 

59.6 

0.3606 

0.5195 

31.61 

30.75 
30.04 
29.45 
28.97 

28.59 
28.28 
28.02 
27.80 

27.59 
27.37 
27.11 
26.79 
26.35 

25.76 


Spot 
Radius 
(10—*  m> 


4.807 

4.258 

3.777 

3.354 

2.981 

2.652 

2.360 

2.100 

1.869 

1.663 

1.478 

1.312 

1.163 

1.028 

0.905 


TABLE  2  (CONTINUED) 

CATHODE  SPOT  PARAMETERS  FOR  A  100  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode  Electric 

Temperature  Field 
(Kelvin)  (10*  V/m) 


3500 

3550 

3600 

3650 

3700 

3750 

3800 

3850 

3900 

3950 


Cathode 

Temperature 

(Kelvin) 


3500 

3550 

3600 

3650 

3700 

3750 

3800 

3850 

3900 

3950 


Electric 
Field 
(10*  V/m) 

E lec  tron 
Current 
Density 
(10*  A/m  =  ) 

Ion 

Current 
Density 
<10“  A/m=> 

Total 
Current 
Density 
(10“  A/m- ) 

3.055 

0.625 

0.753 

1.378 

3.318 

0.841 

0.915 

1.756 

3.598 

1.124 

1.108 

2.232 

3.896 

1.495  I 

1.335 

2.830 

4.214 

1.980 

1.602 

3.582 

4.555 

2.610 

1.917 

4.527 

4.921 

3.430 

2.288 

5.717 

5.315 

4.493 

2.723 

7.216 

5.741 

5.874 

3.236 

9.110 

6.205 

7.670 

3.841 

11.51 

6.712 

10.01 

4.559 

14.57 

7.270 

13.07 

5.412 

IB. 49 

7.889 

17.11 

6.435 

23.54 

8.583 

22.47 

7.670 

30.14 

9.368 

29.69 

9.177 

38.87 

Ion/Total 

Cou  lomb 

Conduc- 

Erosion 

Current 

Density 

Ratio 

Logarithm 

tivity 

(mho/m) 

Rate 

(lO-7- 

atoms/m=-s) 

.5462 

.5212 

.4963 

.4717 


3309 

3489 

3679 

3079 

4091 

4315 

4553 

4000 

5002 

5370 

5701 


3.226 

3.686 

4.193 

4.750 

5.360 

6.025 

6.748 


8.372 

9.276 
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TABLE  3 

CATHODE  SPOT  PARAMETERS  FOR  A  150  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode  Plasma 

Temperature  Temperature 

(Kelvin)  (eV) 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Ca  thode 
T  emperature 
(Kelvin) 


3500 

3550 

3600 

3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Plasma 

Cathode 

T  ransition 

Total  Arc 

T  emperature 

Potential 

Potential 

Potential 

(eV) 

(Volts) 

(Volts) 

(Volts) 

1.48 

22.66 

10.79 

33.45 

1.53 

21.33 

11.41 

32.74 

1.58 

20.15 

12.03 

32.18 

1.64 

19.09 

12.65 

31.74 

1.69 

18.16 

13.25 

31.41 

1.75 

17.33 

13.84 

31.17 

1.81 

16.60 

14.41 

31.01 

1.88 

15.96 

14.93 

30.89 

1.94 

15.41 

15.40 

30.81 

2.02 

14.94 

15.79 

30.73 

2.09 

14.54 

16.08 

30.63 

2.18 

14.22 

16.25 

30.47 

2.27 

13.97 

16.26 

30.23 

2.37 

13.80 

16.06 

29.85 

2.49 

13.70 

15.59 

29.30 

Partial 

Degree  o f 

Return 

Spot 

Pressure 

Ionization 

Coefficient 

Radius 

(atm.) 

<10-*  m) 

10.9 

0.1589 

0.4716 

5.948 

12.6 

0.1674 

0.4742 

5.269 

14.5 

0.1763 

0.4768 

4.674 

16.6 

0.1855 

0.4794 

4.151 

19.0 

0.1952 

0.4821 

3.690 

21.6 

0.2055 

0.4849 

3.282 

24.5 

0.2164 

0.4877 

2.921 

27.6 

0.2281 

0.4907 

2.600 

31.1 

0.2407 

0.4938 

2.315 

34.9 

0.2544 

0.4971 

2.060 

39.1 

0.2696 

0.5006 

1.831 

43.6 

0.2864 

0.5044 

1.626 

48.5 

0.3054 

0.5084 

1.442 

53.8 

0.3271 

0.5129 

1.275 

59.6 

0.3522 

0.5180 

1.123 

'0 
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TABLE  3  (CONTINUED) 

CATHODE  SPOT  PARAMETERS  FOR  A  150  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode  Electric 

Temperature  Field 
(Kelvin)  (10*  Wm> 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Cathode 
T  emperature 
(Kelvin) 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Electric 

Field 

(10*  Wm> 

Electron 
Current 
Density 
(10*  A/m2) 

Ion 

Current 
Density 
(10*  A/m3) 

Total 
Current 
Density 
(10«  A/m3) 

3.011 

0.615 

0.735 

1.349 

3.270 

0.827 

0.893 

1.720 

3.546 

1.105 

1.081 

2.185 

3.840 

1.469 

1.303 

2.771 

4.153 

1.944 

1.564 

3.508 

4.489 

2.562 

1.871 

4.433 

4.850 

3.364 

2.232 

5.596 

5.238 

4.405 

2.657 

7.061 

5.658 

5.755 

3.157 

8.912 

6.115 

7.509 

3.747 

11.26 

6.614 

9.792 

4.445 

14.24 

7.163 

12.78 

5.277 

18.06 

7.772 

16.71 

6.272 

22.98 

8.454 

21.92 

7.473 

29.39 

9.225 

28.92 

8.937 

37.86 

Ion/Total 

Coulomb 

Conduc- 

Erosion 

Current 

Logarithm 

tivity 

Rate 

Density 

Ratio 

(mho/m) 

(lO3'* 

atoms/m3-s) 

0.5443 

10.5 

3279 

3.232 

0.5194 

10.5 

3458 

3.692 

0.4946 

10.5 

3645 

4.200 

0.4700 

10.4 

3842 

4.759 

0.4458 

10.4 

4051 

5.371 

0.4221 

10.4 

4272 

6.038 

0.3989 

10.4 

4506 

6.762 

0.3762 

10.3 

4757 

7.547 

0.3542 

10.3 

5027 

8.392 

0.3329 

10.3 

5318 

9.299 

0.3122 

10.3 

5636 

10.27 

0.2923 

10.3 

5985 

11.30 

0.2729 

10.3 

6372 

12.39 

0.2542 

10.3 

6808 

13.54 

0.2361 

10.3 

7307 

14.74 
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TABLE  4 

CATHODE  SPOT  PARAMETERS  FOR  A  175  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode  Plasma 

Temperature  Temperature 
(Kelvin)  (eV) 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Cathode 

Temperature 

(Kelvin) 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Cathode 

Potential 

(Volts) 


22.61 

21.26 

20.08 

19.04 

18.10 


T  ransition 
Potential 
(Volts) 


Total  Arc 
Potential 
(Volts) 


34.29 

33.63 

33.13 

32.76 


1.75 

17.28 

15.04 

32.32  ' 

1.81 

16.55 

15.66 

32.21 

1.87 

15.92 

16.24 

32.16 

1.94 

15.37 

16.76 

32.13 

2.01 

14.90 

17.21 

32.11 

2.09 

14.50 

17.55 

32.05 

2.17 

14.18 

17.76 

31.95 

2.26 

13.93 

17.80 

31.74 

2.37 

13.76 

17.63 

31.39 

2.48 

13.66 

17.18 

30.84 

Partial 

Degree  of 

Return 

Spot 

Pressure 

Ionization 

Coefficient 

Radius 

(atm.) 

o 

1 

* 

3 

10.9 

0.1578 

0.4714 

6.452 

12.6 

0.1662 

0.4739 

5.711 

14.5 

0.1750 

0.4765 

5 .066 

16.6 

0.1842 

0.4791 

4.499 

19.0 

0.1938 

0.4818 

3.999 

21.6 

0.2040 

0.4845 

3.556 

24.5 

0.2148 

0.4874 

3.166 

27.6 

0.2264 

0.4903 

2.819 

31.1 

0.2389 

0.4934 

2.509 

34.9 

0.2525 

0.4967 

2.233 

39.1 

0.2676 

0.5002 

1.986 

43.6 

0.2843 

0.5039 

1.763 

48.5 

0.3031 

0.5080 

1.563 

53.8 

0.3246 

0.5125 

1.382 

59.6 

0.3495 

0.5175 

1.218 

TABLE  4  (CONTINUED) 

CATHODE  SPOT  PARAMETERS  FOR  A  175  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode 

Temperature 

<Kelvin) 


Electric 
Field 
(10"  V/m) 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Cathode 
T  emperature 
(Kelvin) 


3500 

3550 

3600 

3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


4.468 

4.827 

5.213 

5.631 

6.085 

6.581 

7.128 

7.734 

8.412 

9.178 


Ion/Total 

Current 

Density 

Ratio 


0.5448 

0.5188 

0.4940 

0.4695 

0.4453 

0.4216 

0.3984 

0.3758 

0.3539 

0.3326 

0.3120 

0.2921 

0.2728 

0.2542 

0.2361 


Electron 

Ion 

Total 

Current 

Current 

Current 

Density 

Density 

Density 

<10*  A/m=> 

(10*  A/m=> 

(10*  A/m=> 

0.609 

0.729 

1.338 

0.822 

0.886 

1.708 

1.098 

1.072 

2.171 

1.460 

1.292 

2.752 

1.932 

1.551 

3.483 

2.546 

1.856 

4.402 

3.342 

2.214 

5.557 

4.376 

2.635 

7.011 

5.716 

3.131 

8.847 

7.456 

3.716 

11.17 

9.722 

4.408 

14.13 

12.68 

5.232 

17.91 

16.58 

6.218 

22.79 

21.74 

7.408 

29.15 

28.67 

8.859 

37.53 

Cou  lomb 

Conduc- 

Erosion 

Logarithm 

tivity 

Rate 

(mho/m) 

<10=^ 

atoms/m2-s) 

10.5 

3270 

3.233 

10.5 

3447 

3.694 

10.5 

3634 

4.203 

10.4 

3830 

4.762 

10.4 

4038 

5.374 

10.4 

4258 

6.042 

10.4 

4491 

6.767 

10.3 

4741 

7.552 

10.3 

5009 

8.398 

10.3 

5299 

9.306 

10.3 

5614 

10.28 

10.3 

5961 

11.31 

10.3 

6346 

12.40 

10.3 

6779 

13.55 

10.3 

7274 

14.75 

Figure  5.4  Potentials  for  a  175  Ampere  Arc 
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TABLE  5 

CATHODE  SPOT  PARAMETERS  FOR  A  200  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode 


(Kelvin) 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Cathode 
T  emperature 
(Kelvin) 


Plasma 

Cathode 

T  emperature 

Potential 

(eV) 

(Volts) 

1.48 

22.48 

1.53 

21.21 

1.58 

20.03 

1.63 

18.99 

1.69 

18.06 

1.74 

17.24 

1.80 

16.51 

1.87 

15.88 

1.94 

15.33 

2.01 

14.86 

2.08 

14.47 

2.17 

14.15 

2.26 

13.90 

2.36 

13.73 

2.47 

13.63 

Partial 

Degree  of 

Pressure 

Ionization 

(atm.) 

0.1575 

0.1652 

0.1739 

0.1830 

0.1926 

0.2027 

0.2135 

0.2250 

0.2374 

0.2510 

0.2659 

0.2825 

0.3012 

0.3226 

0.3473 


T  ransition 
Potential 
(Volts) 


12.53 
13.26 
13.99 
14.72 
15.44 
16.15 
16.82 
17.46 
18.03 

18.53 
18.92 
19.17 
19.24 
19.09 
18.65 


Return 

Coefficient 


0.4713 

0.4736 

0.4762 

0.4788 

0.4815 

0.4842 

0.4871 

0.4900 

0.4931 

0.4964 

0.4998 

0.5036 

0.5076 

0.5121 

0.5171 


Total  Arc 
Potential 
(Volts) 


35.01 

34.47 

34.02 

33.71 

33.50 

33.38 
33.34 
33.34 
33.37 

33.39 

33.39 
33.32 
33.15 
32.82 
32.28 


Spot 
Radius 
(10—  m) 


5. 
4.823 


3.395 

3.022 

2.690 

2.394 
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TABLE  5  (CONTINUED) 

CATHODE  SPOT  PARAMETERS  FOR  A  200  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode 

Temperature 

(Kelvin) 


3500 

3550 

3600 

3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 

Cathode 
Tempera ture 
(Kelvin) 


3500 

3550 

3600 

3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Electric 

Field 

<10*  0/m) 


2.990 

3.241 

3.515 

3.806 

4.117 

4.450 

4.808 

5.193 

5.609 

6.061 

6.556 

7.100 

7.703 

8.378 

9.141 

Ion/Total 

Current 

Density 

Ratio 

0.5438 

0.5183 

0.4935 

0.4690 

0.4449 

0.4212 

0.3981 

0.3755 

0.3536 

0.3323 

0.3118 

0.2919 

0.2727 

0.2541 

0.2360 


Electron 
Current 
Density 
(10*  A/m=> 

0.610 

0.818 

1.093 

1.453 

1.923 

2.533 

3.326 

4.353 

5.685 

7.414 

9.665 

12.61 

16.47 
21.60 

28.47 

Coulomb 

Logarithm 


Ion 

Current 

Density 
(10*  A/m=) 

0.727 

0.880 

1.065 

1.284 

1.541 

1.844 

2.199 

2.617 

3.110 

3.691 

4.379 

5.197 

6.176 

7.357 

8.797 

Conduc¬ 

tivity 

(mho/m) 


3267 

3439 

3625 

3820 

4027 

4246 

4479 

4727 

4994 

5283 

5597 

5942 

6325 

6755 

7247 


Total 
Current 
Density 
(10*  A/m2) 

1.338 

I. 698 
2.159 
2.737 
3.464 

4.377 
5.525 
6.970 
8.795 

II. 11 
14.04 
17.80 
22.65 
28.95 

37.27 

Erosion 

Rate 

(10:^,,• 

atoms/m3-s) 

3.234 

3.696 

4.205 

4.765 

5.377 
6.045 
6.771 
7.556 
8.403 
9.312 

10.28 
11.32 
12.41 
13.56 
14.77 
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TABLE  6 

CATHODE  SPOT  PARAMETERS  FOR  A  300  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode 

Plasma 

Cathode 

T  ransition 

Total  Arc 

Temperature 

Temperature 

Potential 

Potential 

Potential 

(Kelvin) 

(eV) 

(Volts) 

(Volts) 

(Volts) 

3500 

1.47 

22.32 

15.45 

37.77 

3550 

1.52 

21.06 

16.36 

37.42 

3600 

1.57 

19.90 

17.28 

37.17 

3650 

1.62 

18.86 

18.20 

37.06 

3700 

1.68 

17.94 

19.11 

37.05 

3750 

1.74 

17.13 

20.00 

37.13 

3800 

1.80 

16.41 

20.87 

37.28 

3850 

1.86 

15.78 

21.69 

37.47 

3900 

1.93 

15.24 

22.45 

37.69 

3950 

2.00 

14.77 

23.11 

37.89 

4000 

2.07 

14.38 

23.66 

38.04 

4050 

2.15 

14.07 

24.04 

38.11 

4100 

2.25 

13.82 

24.23 

38.05 

4150 

2.34 

13.65 

24.15 

37.80 

4200 

2.46 

13.55 

23.75 

37.30 

Cathode 

Partial 

Degree  of 

Return 

Spot 

T  emperature 

Pressure 

Ionization 

Coefficient 

Radius 

(Kelvin) 

(atm.) 

(10—*  m) 

3500 

10.9 

0.1550 

0.4706 

8.513 

3550 

12.6 

0.1626 

0.4730 

7.554 

3600 

14.5 

0.1712 

0.4755 

6.701 

3650 

16.6 

0.1801 

0.4781 

5.951 

3700 

19.0 

0.1896 

0.4808 

5.290 

3750 

21.6 

0.1995 

0.4835 

4.707 

3800 

24.5 

0.2101 

0.4863 

4.189 

3850 

27.6 

0.2214 

0.4892 

3.730 

3900 

31.1 

0.2336 

0.4923 

3.321 

3950 

34.9 

0.2470 

0.4955 

2.956 

4000 

39.1 

0.2616 

0.4990 

2.629 

4050 

43.6 

0.2779 

0.5027 

2.336 

4100 

48.5 

0.2963 

0.5067 

2.071 

4150 

53.8 

0.3172 

0.5111 

1.832 

4200 

59.6 

0.3415 

0.5160 

1.616 

TABLE  6  (CONTINUED) 

CATHODE  SPOT  PARAMETERS  FOR  A  300  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode 

Temperature 

(Kelvin) 

Electric 
Field 
(10*  V/m> 

Elec  tron 
Current 
Density 
<10a  A/m=> 

Ion 

Current 
Density 
(10*  A/m=> 

Total 
Current 
Density 
<10*  A/m=> 

3500 

2.959 

0.603 

0.715 

1.318 

3550 

3.207 

0.808 

0.865 

1.673 

3600 

3.478 

1.080 

1.047 

2.127 

3650 

3.767 

1.435 

1.261 

2.696 

3700 

4.075 

1.898 

1.514 

3.412 

3750 

4.404 

2.500 

1.811 

4.311 

3800 

4.758 

3.280 

2.161 

5.441 

3850 

5.139 

4.292 

2.572 

6.863 

3900 

5.551 

5.603 

3.055 

9 .658 

3950 

5.998 

7.304 

3.625 

10.93 

4000 

6.487 

9.516 

4.300 

13.82 

4050 

7.025 

12.40 

5.103 

17.51 

4100 

7.621 

16.20 

6.063 

22.26 

4150 

8.287 

21.22 

7.221 

28.44 

4200 

9.040 

27.95 

8.632 

36.58 

Cathode 

Ion/Total 

Coulomb 

Conduc- 

Erosion 

Temperature 

(Kelvin) 

Current 

Density 

Ratio 

Logarithm 

tivity 

(mho/m) 

Rate 

(1027, 

atoms/m2-s) 

3500 

0.5423 

10.5 

3246 

3.238 

3550 

0.5169 

10.5 

3416 

3.701 

3600 

0.4922 

10.5 

3600 

4.211 

3650 

0.4678 

10.4 

3794 

4.771 

3700 

0.4438 

10.4 

3999 

5.385 

3750 

0.4202 

10.4 

4216 

6.054 

3800 

0.3971 

10.4 

4446 

6.781 

3850 

0.3747 

10.3 

4692 

7.568 

3900 

0.3528 

10.3 

4956 

8.417 

3950 

0.3317 

10.3 

5241 

9.328 

4000 

0.3113 

10.3 

5551 

10.30 

4050 

0.2915 

10.3 

5891 

11.34 

4100 

0.2724 

10.3 

6269 

12.44 

4150 

0.2539 

10.3 

6693 

13.59 

4200 

0.2360 

10.3 

7176 

14.80 

TABLE  7 

CATHODE  SPOT  PARAMETERS  FOR  A  400  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode 

Temperature 

(Kelvin) 

Plasma 

T  emperature 

<eV) 

Cathode 

Potential 

(Volts) 

T  ransition 
Potential 
(Volts) 

Total  Arc 
Potential 
(Volts) 

3500 

1.47 

22.22 

17.91 

40.13 

3550 

1.52 

20.97 

18.97 

39.94 

3600 

1.57 

19.81 

20.05 

39.86 

3650 

1.62 

18.78 

21.12 

39.91 

3700 

1.67 

17.87 

22.20 

40.07 

3750 

1.73 

17.06 

23.25 

40.31 

3800 

1.79 

16.35 

24.28 

40.63 

3850 

1.85 

15.72 

25.25 

40.98 

3900 

1.92 

15.18 

26.16 

41.34 

3950 

1.99 

14.72 

26.97 

41.69 

4000 

2.07 

14.33 

27.64 

41.98 

4050 

2.15 

14.02 

28.15 

42.16 

4100 

2.24 

13.77 

28.42 

42.20 

4150 

2.34 

13.60 

28.41 

42.01 

4200 

2.45 

13.50 

28.03 

41.53 

Cathode 

Temperature 

(Kelvin) 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Partial 

Pressure 

(atm.) 


Degree  oi 
Ionization 


0.1535 

0.1610 

0.1695 

0.1784 

0.1877 

0.1976 

0.2081 

0.2193 

0.2314 

0.2446 

0.2591 

0.2752 

0.2934 

0.3141 

0.3380 


Return 

Coefficient 


0.4702 

0.4726 

0.4751 

0.4777 

0.4803 

0.4830 

0.4858 

0.4888 

0.4918 

0.4950 

0.4984 

0.5021 

0.5061 

0.5105 

0.5154 


Spot 
Radius 
(10-*  m> 


9.874 

8.762 

7.772 

6.903 

6.136 

5.459 

4.860 

4.327 

3.853 

3.430 

3.051 

2.710 

2.404 

2.127 

1.876 
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TABLE  7  (CONTINUED) 

CATHODE  SPOT  PARAMETERS  FOR  A  400  AMPERE  ARC  ON  A 

COPPER  CATHODE 


Cathode 

Temperature 

(Kelvin) 


3650 

3700 

3750 

3800 

3850 

3900 

3950 

4000 

4050 

4100 

4150 

4200 


Cathode 

Temperature 

(Kelvin) 


Electric 

Field 

(10«  V/m) 

Electron 
Current 
Density 
(10*  A/m2) 

Ion 

Current 
Density 
(10«  A/m2) 

2.940 

0.599 

0.707 

3.187 

0.803 

0.856 

3.456 

1.072 

1.036 

3.743 

1.424 

1.248 

4.049 

1.883 

1.498 

4.377 

2.480 

1.792 

4.728 

3.253 

2.138 

5.107 

4.256 

2.544 

5.516 

5.554 

3.022 

5.961 

7.238 

3.586 

6.446 

9.427 

4.254 

6.980 

12.28 

5.047 

7.572 

16.03 

5.996 

8.234 

21.00 

7.141 

8.981 

27.64 

8.534 

Ion/Total 

Coulomb 

Conduc- 

Current 

Density 

Ratio 

Logarithm 

tivity 

(mho/m) 

Total 
Current 
Density 
(10«  A/m=> 


2.672 

3.3822 

4.272 

5.391 

6.800 

8.577 

10.83 

13.68 

17.33 

22.03 

28.14 

36.17 


Erosion 

Rate 

(1027 

atoms/m2-s) 


3.240 
3.704 
4.214 
4.775 
5.389 
6.059 
6.787 
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TABLE  8 

CATHODE  SPOT  PARAMETERS  FOR  STABLE  ARCS 


Arc  Current 


Parameter 

200 

300 

400 

Cathode  Temperature 
(Kelvin) 

3850 

3700 

3600 

Plasma  Temperature  (eV) 

L87 

L68 

1.57 

Cathode  Potential  (Volts) 

15.88 

17.94 

19.81 

Spot  Radius  (10~*  m) 

3.022 

5.290 

7.772 

Ion  Current  Density  (108 

A/m2 ) 

2.617 

1.514 

L036 

Electron  Current  Density 
(10»  A/m2) 

4.353 

1.898 

L072 

Total  Current  Density  (10 8 
A/m2 ) 

6.970 

3.412 

2.108 

Erosion  Rate  (10 2  v 
atoms /m2 -s) 

7.556 

5.385 

4.214 

’•I 
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This  study  has  attempted  to  answer  some  of  the  questions 
associated  with  the  cathode  phenomena  in  a  unipolar  arc. 
However,  a  more  detailed  analysis  of  the  cathode  spot  must  be 
accomplished  before  a  consistent  understanding  of  the  diverse 
phenomena  occurring  within  the  arc  can  be  achieved.  The 
stationary  model  presented  her  is  only  a  beginning  and  a  tool 
for  the  development  of  future  models.  It  can  not  answer  the 
questions  regarding  arc  initiation,  extinction,  and  spot  motion 
that  will  have  to  be  addressed  in  future  models  if  they  are  to 
be  accepted.  Additional  research  with  stationary  models  should 
consider  including  the  reverse  electron  current  from  the  plasma 
back  to  the  metal,  changing  from  a  constant  cathode  temperature 
in  the  arc  radius  to  a  constant  heat  flux  with  a  varying 
temperature,  considering  the  ohmic  heating  of  the  cathode  spot 
by  the  electron  current,  and  varying  the  cathode  material 
parameters  to  determine  their  influence  on  the  arc  stability.  It 
is  hoped  that  this  study  will  help  to  stimulate  further  interest 
and  research  into  the  cathode  spot  of  a  unipolar  arc. 
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APPENDIX  A  GENERAL  EMISSION  PROGRAM 


4c  *  *  **  *  *  *  **  *  *****  *  ******  **  *  **  *  *  ****  *  ****  *  ********  * 
* 

EMM1.EXE 


PURPOSE 

CALCULATE  THE  ELECTRON  CURRENT  DENSITY 
EMITTED  FOR  A  RANGE  OF  ELECTRIC  FIELDS  AT 
THE  SURFACE  OF  A  CATHODE  WITH  A  GIVEN  WORK 
FUNCTION  AND  CATHODE  TEMPERATURE. 


DESCRIPTION  OF  VARIABLES 


c 

* 

T 

CATHODE  TEMPERATURE 

c 

* 

F 

ELECTRIC  FIELD 

c 

* 

FO 

INITIAL  ELECTRIC  FIELD 

c 

* 

PHI 

WORK  FUNCTION 

c 

* 

DF 

ELECTRIC  FIELD  INCREMENT 

c 

* 

J 

CURRENT  DENSITY 

c 

* 

ECODE 

ERROR  CODE 

c 

* 

0  =  GOOD 

c 

* 

1  =  ERROR 

c 

* 

TYP 

TYPE  OF  EMISSION 

c 

* 

1=  THERMIONIC 

c 

* 

2=  INTERMEDIATE 

c 

* 

3=  FIELD 

c 

* 

4=  NONE 

c 

* 

c 

* 

SUBROUTINES 

REQUIRED 

c 

* 

EMMAPP(T,  F,  PHI.TYP) 

c 

* 

TFEMMKT,  F,  PHI,  J,  ECODE) 

c 

* 

IFTEMM1(T,F, PHI,  J, ECODE) 

c 

* 

FTEMMKT,  F,  PHI,  J,  ECODE) 

************************************************* 


REAL*8  T,F,FO,PHI,DF,J 
INTEGER  I,N,ECODE,TYP 
OPEN(4 ,  FILE= 'CON') 

OPEN(6  ,FILE='EMML  DAT',  ST  ATUS  =  'NEW') 

WRITE**,*)  'ENTER  THE  NUMBER  OF  DATA  POINTS:  ' 
READ**,*)  N 

WRITE**,*)  'ENTER  THE  ABSOLUTE  TEMPERATURE  (K): 
READ**,*)  T 

WRITE**,*)  'ENTER  INITIAL  ELECTRIC  FIELD  (V/M):  ' 
READ**,*)  FO 

WRITE**,*)  'ENTER  ELECTRIC  FIELD  INCREMENT  (V/M): 
READ**,*)  DF 
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WRITE(*,*)  'ENTER  WORK  FUNCTION  (eV):  ' 

READ(*,*)  PHI 

F=F0 

J=0. 

DO  20  1=1, N 

CALL  EMMAPP(T,  F,  PHI,  T YP) 

IF(TYP.EQ.l)  CALL  TFEMM1(T,F,PHI, J.ECODE) 
IFfTYP.EQ.2)  CALL  IFTEMM1(T,F,  PHI,  J.ECODE) 
IF(TYP.EQ.3)  CALL  FTEMMKT.F, PHI, J.ECODE) 
IF(TYP.EQ.4)  GOTO  10 
WRITE(6,900)  F,  J,  ECODE,  TYP 
F=F  +  DF 
CONTINUE 
STOP 

FORMATC1X,  D12.6,5X,D12.6,5X,I5,5X,I5) 

END 
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APPENDIX  B  GENERAL  EMISSION  SOBROOTINB 

************************************************* 

* 

*  SUBROUTINE  EMM 

* 

*  PURPOSE 

*  CALCULATE  THE  ELECTRON  CURRENT  DENSITY 

*  FOR  A  GIVEN  ELECTRIC  FIELD  AT  THE  SURFACE 

*  OF  A  CATHODE  WITH  A  GIVEN  WORK  FUNCTION  AND 

*  CATHODE  TEMPERATURE. 

* 

*  DESCRIPTION  OF  VARIABLES 

*  T  CATHODE  TEMPERATURE 

*  F  ELECTRIC  FIELD 

*  PHI  WORK  FUNCTION 

*  J  CURRENT  DENSITY 

*  ECODE  ERROR  CODE 

*  0  =  GOOD 

*  1  =  ERROR 

*  TYP  TYPE  OF  EMISSION 

*  1=  THERMIONIC 

*  2=  INTERMEDIATE 

*  3  —  FIELD 

*  4=  NONE 

* 

*  SUBROUTINES  REQUIRED 

*  EMMAPP(T,  F,  PHI,  TYP) 

*  TFEMM1(T,F,PHI,J,  ECODE) 

*  IFTEMM1(T,F,  PHI,  J,  ECODE) 

*  FTEMMKT.F,  PHI,  J,  ECODE) 

* 

************************************************* 

SUBROUTINE  EMM(T,F,PHI, J, ECODE, TYP) 

REAL*8  T.F.PHI.J 
INTEGER  ECODE, TYP 
CALL  EMMAPP(T,  F,  PHI,  T YP) 

IF(TYP.EQ.  1)  CALL  TFEMM1(T, F, PHI, J, ECODE) 

IF(TYP.EQ.2)  CALL  IFTEMMKT, F,  PHI, J, ECODE) 

IF(TYP.EQ.3)  CALL  FTEMM1(T, F, PHI, J, ECODE) 

RETURN 

END 
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APPIHDIZ  C  miSSIOM  APPLICABILITY  SOBBOOTIMK 

************************************************** 

* 

*  SUBROUTINE  EMMAPP 

* 

*  PURPOSE 

*  DETERMINES  APPLICABILITY  OF  THE  SUBROUTINES 

*  TFEMM1,  IFTEMM1,  AND  FTEMM1 


DESCRIPTION 

T 

FO 

PHIO 

TYP 


SUBROUTINES 

NONE 


OF  VARIABLES 
CATHODE  TEMPERATURE 
ELECTRIC  FIELD  AT  CATHODE  SURFACE 
WORK  FUNCTION 
TYPE  OF  EMISSION 
1=  THERMIONIC 
2=  INTERMEDIATE 
3=  FIELD 
4=  NONE 

REQUIRED 


SUBROUTINE  EMMAPP(T,FO, PHIO, TYP) 

RE AL*8  T, F, PHI, M, E, K, H, EPS, PI, PHIO , D, D1 

RE AL*8  FO , El, FI, KTO , KT, Y, YE,TY, TE, V, VE, ETA, ET AO 

REAL*8  TEST1,TEST2,TEST3,TEST4 

INTEGER  ECODE1,  ECODE2 ,  ECODE3 ,  T YP,  N 

M=9.1095E-31 

E=1.6022E-19 

K=L3807E-23 

H=L0546E-34 

EPS  =  8.8542E-1 2 

PI=3. 141592654 

E1=(4*PI*EPS*H)**2/(M*E**4) 

F1=E1**2*E**3/(4*PI*EPS) 

F=F0*F1 

KT0=K*T 

KT=KT0*E1 

PHI=PHI0*E1*E 

N=0 

D=F**(3./4.)/(PI*KT) 

IF(l-D)  20,20,10 

TEST1=  DLOG((l-  D)/D)-  1/(D*(1-  D)) 

TEST2  =  -PI/(F**(3./4.))*(PHI-  DSQRT(F)) 

TEST3  =  DLOG((l-  D)/D)- 1/(1-  D) 
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TEST4=-PI/(F**(l./8.)) 

IF((TESTi.GT.TEST2).AND.(TEST3.GT.TEST4))  THEN 
ECODE1=0 
ELSE 

EC0DE1=1 
ENDIF 
GOTO  30 
EC0DE1=1 

ETA0  =  -F**2/(8*KT**2) 

Y=DSQRT(F)/PHI 
CALL  VTY(Y,  V,TY,  IER1) 

YE=DSQRT(F)/(-ETAO) 

CALL  VTY(YE,VE,TE,IER2) 
ETA=-F**2/(8*KT**2*TE**2) 

IF(DABS(ETA-ETAO).GT.  LOD-6)  THEN 
ETA0  =  ETA-(ETA-ETA0)/2.1 
N=N+1 

IF(N.GT.200)  THEN 
ECODE2=l 
GOTO  50 
ENDIF 
GOTO  40 
ENDIF 

D=2*DSQRT(2/YE)*TE/PI 
TEST1=  1/YE 

TEST2  =  1  +  F**(L/4.)*D/(PI*(D-1)) 

TEST3  =  -F**2/(8*KT**2*TE**2) 

TEST4  =  -PHI+KT/(1-F/(2*DSQRT(2*PHD*KT*TY)) 
IF((TEST1.GT.TEST2).AND.(TEST3.GT.TEST4))  THEN 
ECODE2=0 
ELSE 

ECODE2=l 

ENDIF 

D=  (2  *  DS  QRT(2 . 0  *  PHI)*  PI*  KT  )/F 
D1=4*DSQRT(2*PHI)*PHI/(3*F) 

TEST1=PHI-  DSQRT(F) 

TEST2  =  F**(3./4.)/PI  +  KT/(1-D*TY/PI) 

TEST3  =  1-D*TY/PI 

TEST4=KT*DSQRT((V*PHI/F)*DSQRT(2*PHI)/(PHI**2-F)) 
IF((TESTLGT.TEST2).AND.(TEST3.GT.TEST4))  THEN 
ECODE3=0 
ELSE 

ECODE3=l 

ENDIF 

IF(ECODELEQ.O)  then 
TYP=1 
ELSE 

IF(ECODE2.EQ.O)  THEN 
TYP=2 
ELSE 

IF(ECODE3.EQ.O)  THEN 
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TYP=3 

ELSE 

TYP=4 

ENDIF 

ENDIF 

ENDIF 

RETURN 

END 
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PURPOSE 

CALCULATES  THE  ELECTRON  CURRENT  DENSITY 
FOR  THERMIONIC  FIELD  EMISSION 

DESCRIPTION  OF  VARIABLES 


APPENDIX  D  EMISSION  SOBBOOTIMR 

************************************************** 

* 

*  SUBROUTINE  TFEMM1 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

************************************************** 


T 

FO 

PHIO 

JO 

ECODE 


SUBROUTINES 

NONE 


CATHODE  TEMPERATURE 
ELECTRIC  FIELD  AT  CATHODE  SURFACE 
WORK  FUNCTION 
CURRENT  DENSITY 
ERROR  CODE 
0=  GOOD 
1=  ERROR 

REQUIRED 


10 


S UBROUTINE  TFEMM1(T, FO , PHIO , JO , ECODE) 

REAL*8  T, F, PHI, J, M, E, K, H, EPS, PI, PHIO , D 

REAL*8  FO , JO , El, FI, Jl, KTO , KT 

REAL*8  TEST1,TEST2,TEST3,TEST4 

INTEGER  ECODE 

M-9.1095E-31 

E=L6022E-19 

K=L3807E-23 

H=L0546E-34 

EPS=8.8542E-12 

PI=3. 141592654 

E1=(4*PI*EPS*H)**2/(M*E**4) 

F1=E1**2*E**3/(4*PI*EPS) 

J1=E1**2*H**3/(M*E) 

F=F0*F1 

KT0=K*T 

KT=KT0*E1 

PHI=PHI0*E1*E 

D=F**(3./4.)/(PI*KT) 

J=(KT**2/(2*PI**2))*((PI*D)/DSIN(PI*D))* 
1  DEXP(-(PHI-  DSQRT(F))/KT) 

J0  =  J/J1 

IF(l-D)  20,20,10 
TEST1=  DLOG((l-  D)/D)-  1/(D*(1-  D» 
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TEST2= -  PI/(F**(3.  /4.  ))*(PHI-  DSQRT(F)) 
TEST3=DL0G((1-  D)/D)- 1/(1-  D) 

TEST4 = -  PI/(F**(L  /8 . )) 

IF((TESTLGT.TEST2).AND.(TEST3.GT.TEST4))  THEN 
ECODE=0 
ELSE 

EC0DE=1 

ENDIF 

RETURN 

EC0DE=1 

RETURN 

END 
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APPENDIX  E  g  THLD- THERMION IC  EMISSION  SOBROOTINE  < 


************************************************** 

* 

*  SUBROUTINE  FTEMM1 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


PURPOSE 

CALCULATES  THE  ELECTRON  CURRENT  DENSITY 
FOR  FIELD-THERMIONIC  EMISSION 

DESCRIPTION  OF  VARIABLES 

T  CATHODE  TEMPERATURE 

FO  ELECTRIC  FIELD  AT  CATHODE  SURFACE 

PHIO  WORK  FUNCTION 

JO  CURRENT  DENSITY 

ECODE  ERROR  CODE 

0=  GOOD 
1=  ERROR 

SUBROUTINES  REQUIRED 
VTY(Y,  V.TY.IER) 


************************************************** 


SUBROUTINE  FTEMM1(T,F0, PHIO,  JO, ECODE) 

REAL*8  T, F, PHI, J,  M, E, K, H, EPS, PI, PHI 0 , D, Dl, T Y 

REAL*8  TEST1,TEST2,TEST3,TEST4, Y, V,C1,C2,C3 

REAL*8  FO ,  JO ,  El,  FI,  Jl,  KTO , KT 

INTEGER  ECODE 

M=9.1095E-31 

E=1.6022E-19 

K=1.3807E-23 

H=  1.0546E-34 

EPS=8.8542E-12 

PI=3. 141592654 

E1=(4*PI*EPS*H)**2/(M*E**4) 

F1=E1**2*E**3/(4*PI*EPS) 

J1=E1**2*H**3/(M*E) 

F=F0*F1 

KT0=K*T 

KT=KT0*E1 

PHI=PHI0*E1*E 

D=(2*DSQRT(2.0*  PHI)*  PI*  KT  )/F 
D1=4*DSQRT(2*PHI)*PHI/(3*F) 

Y=DSQRT(F)/PHI 
CALL  VTY(Y  V.TY.IER) 

J=F**2*(D*Ty'/DS  IN(D*TY))*DEXP(-D1*V)/(16*PI**2*PHI*TY**2) 
J0  =  J/J1 


TEST1=  PHI-  DSQRT(F) 

TEST2=F**(3./4.)/PI  +  KT/(1-D*TY/PI) 

TEST3=1-  D*TY/PI 

TEST4=KT*DSQRT((V*PHI/F)*DSQRT(2*PHD/(PHI**2-F)) 
IF((TESTLGT.TEST2).AND.(TEST3.GT.TEST4))  THEN 
ECODE=0 
ELSE 

EC0DE=1 

ENDIF 

RETURN 

END 
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APPENDIX  F  COMPLETE  ELLIPTIC  INTEGRAL  Ql_  XflE.  1ST 

EIMH 


************************************************** 
SUBROUTINE  CEL1 
PURPOSE 

CALCULATE  COMPLETE  ELLIPTIC  INTEGRAL  OF  FIRST 
KIND 

USAGE 

CALL  CEL1(RES,  AK.IER) 

DESCRIPTION  OF  PARAMETERS 
RES  -  RESULT  VALUE 

AK  -  MODULUS  (INPUT) 

IER  -  RESULTANT  ERROR  CODE  WHERE 

IER=0  NO  ERROR 

IER=1  AK  NOT  IN  RANGE  -1  TO  +1 

REMARKS 

THE  RESULT  IS  SET  TO  LE75  IF  ABS(AK)  GE  1 
FOR  MODULUS  AK  AND  COMPLEMENTARY  MODULUS  CK, 
EQUATION  AK*AK+CK*CK=LO  IS  USED. 

AK  MUST  BE  IN  THE  RANGE  -1  TO  +1 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

METHOD 

DEFINITION 

CEL1(AK)=INTEGRAL(1/SQRT((1+T*T)*(1+(CK*T)**2)), 

SUMMED  OVER  T  FROM  0  TO  INFINITY). 

EQUIVALENT  ARE  THE  DEFINITIONS 
CELl(AK)=INTEGRAL(l/(COS(T)SQRT(l+(CK*TAN(T))**2)) 
SUMMED  OVER  T  FROM  0  TO  PI/2), 
CEL1(AK)=INTEGRAL(1/SQRT(1-(AK*SIN(T))**2), 

SUMMED  OVER  T  FROM  0  TO  PI/2),  WHERE 
K=SQRT(L-CK*CK). 

EVALUATION 

LANDENS  TRANSFORMATION  IS  USED  FOR  CALCULATION. 
REFERENCE 

R.BULIRSCH,  'NUMERICAL  CALCULATION  OF  ELLIPTIC 
INTEGRALS  AND  ELLIPTIC  FUNCTIONS',  HANDBOOK  SERIES 
SPECIAL  FUNCTIONS,  NUMERISCHE  MATHEMATIK  VOL.  7, 
1965,  PP.  78-90. 
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************************************************** 

SUBROUTINE  CELKRES, AK.IER) 

REAL*8  ARI,GEO,AK,RES, AARI 

IER=0 

ARI=2. 

GEO=(0.5-AK)+0.5 

GEO=GEO+GEO*AK 

RES=0.5 

IF(GEO)l,2,4 

IER=1 

RES=LE38 

RETURN 

GEO=GEO*AARI 

GEO=DSQRT(GEO) 

GEO=GEO+GEO 

AARI=ARI 

ARI=ARI+GEO 

RES=RES+RES 

IF(GEO/AARI-0. 9999)3, 5,5 

RES=RES/ARI*6.283185E0 

RETURN 

END 
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APPENDIX  6  COMPLETE  ELLIPTIC  INTEGRAL  OF  2ND  EINn 


************************************************** 
SUBROUTINE  CEL2 
PURPOSE 

COMPUTES  THE  GENERALIZED  COMPLETE  ELLIPTIC 
INTEGRAL  OF  SECOND  KIND. 

USAGE 

CALL  CEL2(RES,AK,  A.B.IER) 

DESCRIPTION  OF  PARAMETERS 
RES  -  RESULT  VALUE 

AK  -  MODULUS  (INPUT) 

A  -  CONSTANT  TERM  IN  NUMERATOR 

B  -  FACTOR  OF  QUADRATIC  TERM  IN  NUMERATOR 

IER  -  RESULTANT  ERROR  CODE  WHERE 

IER=0  NO  ERROR 

IER=1  AK  NOT  IN  RANGE  -1  TO  +1 

REMARKS 

FOR  ABS(AK)  GE  1  THE  RESULT  IS  SET  TO  LE75  IF  B 
IS  POSITIVE,  TO  -LE75  IF  B  IS  NEGATIVE. 

SPECIAL  CASES  ARE 

K(K)  OBTAINED  WITH  A  =  1,  B  =  1 

E(K)  OBTAINED  WITH  A  =  1,  B  =  CK*CK  WHERE  CK  IS 

COMPLEMENTARY  MODULUS. 

B(K)  OBTAINED  WITH  A  =  1,  B  =  0 

D(K)  OBTAINED  WITH  A  =  0,  B  =  1 

WHERE  K,  E,  B,  D  DEFINE  SPECIAL  CASES  OF  THE 

GENERALIZED  COMPLETE  ELLIPTIC  INTEGRAL  OF  SECOND 

KIND  IN  THE  USUAL  NOTATION,  AND  THE  ARGUMENT  K  OF 

THESE  FUNCTIONS  MEANS  THE  MODULUS. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

METHOD 

DEFINITION 

RES  =  INTEGRAL((A+B*T*T)/(SQRT((1+T*T) 
*(1+(CK*T)**2))*(1+T*T)) 

SUMMED  OVER  T  FROM  0  TO  INFINITY). 

EVALUATION 

LANDENS  TRANSFORMATION  IS  USED  FOR  CALCULATION. 
REFERENCE 


noooooo 


t,  VuVOTwrawr : 


R.BULIRSCH,  'NUMERICAL  CALCULATION  OF  ELLIPTIC 
INTEGRALS  AND  ELLIPTIC  FUNCTIONS',  HANDBOOK  SERIES 
SPECIAL  FUNCTIONS,  NUMERISCHE  MATHEMATIK  VOL.  7, 
1965,  PP.  78-90. 

*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  **  ***  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ** 

SUBROUTINE  CEL2(RES,AK, A,B,IER) 

REAL*8  RES, AK, A, B, ARI, GEO, Al, BO , AARI 

IER=0 

ARI=2. 

GEO= (0 . 5  -  AK)+0 . 5 

GEO=GEO+GEO*AK 

RES=A 

A1=A+B 

B0  =  B+B 

IF(GEO)l,2,6 

1  IER=1 

2  IF(B)3,8,4 

3  RES=-LE38 
RETURN 

4  RES=LE38 
RETURN 

5  GEO=GEO*AARI 

6  GEO=  DSQRT(GEO) 

GEO=GEO+GEO 

AARI=ARI 

ARI=ARI+GEO 

BO  =  BO +RES*GEO 
RES=A1 
B0  =  B0+B0 
A1=B0/ARI+A1 
IF(GEO/AARI- 0.9999)5, 7, 7 

7  RES=A1/ARI 
RES=RES+0.5707963E0*RES 

8  RETURN 
END 
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APPKMDIX  H 


ELLIPTIC  gDMCTIOHS 


************************************************** 

* 

*  SUBROUTINE  VTY 

* 

*  PURPOSE 

*  CALCULATES  THE  FOWLER- NORDHIEM  ELLIPTIC 

*  FUNCTIONS 

* 

*  DESCRIPTION  OF  VARIABLES 

*  Y  FUNCTION  ARGUMENT 

*  V  FOWLER- NORDHIEM  ELLIPTIC  V 

*  TY  FOWLER-NORDHIEM  ELLIPTIC  T 

*  IER  ERROR  CODE 

*  0=  GOOD 

*  1=  ERROR 

* 

*  SUBROUTINES  REQUIRED 

*  CEL1(RES,AK,IER) 

*  CEL2(RES,AK,A,  B,IER) 

* 

************************************************** 

SUBROUTINE  VTY(Y, V,TY,IER) 

REAL*8  Y, V,TY,C1,C2,C3 
INTEGER  IER.IER1.IER2 
IF(Y.GT.  1.0)  THEN 

Cl=  DSQRT((Y- 1)/(2  *  Y)) 

CALL  CEL1(C2,C1, IER1) 

CALL  CEL2(C3,C1,L0D00,1.0D00-C1*C1,IER2) 
V=-DSQRT(Y/2)*(-2*C3  +  (Y+1)*C2) 
TY=1/DSQRT(1+Y)*((1+Y)*C3  -  Y*C2) 

ELSE 

Cl=  DSQRT((1- Y)/(l+Y)) 

CALL  CEL1(C2,C1,IER1) 

CALL  CEL2(C3 , Cl,  1.0  DO 0 , L 0 DO 0 -  C1*C1, IER2 ) 

V=  DSQRT(1+  Y)*(C3  - Y*C2 ) 

TY=1/DSQRT(1+Y)*((1+Y)*C3  -Y*C2) 

ENDIF 

IF((IERLEQ.l).OR.(IER2.EQ.l))  THEN 
IER=1 
ELSE 
IER=0 
ENDIF 
RETURN 
END 


APPENDIX  I 


=J_  EMISSION 


************************************************** 

* 

*  SUBROUTINE  IFTEMM1 


PURPOSE 

CALCULATES  THE  ELECTRON  CURRENT 
FOR  INTERMEDIATE  FIELD- THERMIONIC 

DESCRIPTION  OF  VARIABLES 


DENSITY 

EMISSION 


SURFACE 


T  CATHODE  TEMPERATURE 

FO  ELECTRIC  FIELD  AT  CATHODE 

PHIO  WORK  FUNCTION 

JO  CURRENT  DENSITY 

ECODE  ERROR  CODE 

0=  GOOD 
1=  ERROR 


*  SUBROUTINES  REQUIRED 

*  VTY(Y,  V,TY,IER) 

* 

************************************************** 

SUBROUTINE  IFTEMM1(T, FO , PHIO , JO , ECODE) 

REAL*8  T,  F,  PHI, J,  M,  E,  K,  H,  EPS ,  PI,  PHIO , D, TE, T Y 

REAL*8  TEST1,TEST2,TEST3,TEST4, Y, V,C1,C2,C3, VE 

RE AL*8  FO , JO , El, FI, Jl, KTO , KT, ET A, ETAO , YE, THETA 

INTEGER  ECODE 

M=9.1095E-31 

E=1.6022E-19 

K=L3807E-23 

H=L0546E-34 

EPS=8.8542E-12 

PI=3. 141592654 

E1=(4*PI*EPS*H)**2/(M*E**4) 

F1=E1**2*E**3  /(4*PI*EPS) 

Jl=El**2*H**3  /(M*E) 

F=F0*F1 

KT0=K*T 

KT=KT0*E1 

PHI=PHI0*E1*E 

ETA0  =  -F**2/(8*KT**2) 

Y=DSQRT(F)/PHI 
CALL  VTY(Y,V,TY,IER1) 

YE=  DSQRT(F)/(-  ETAO ) 

CALL  VTY(YE,VE,TE,IER2) 

ETA=-F**2/(8*KT**2*TE**2) 
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IF(DABS(ETA-ETAO ).GT.  LOD-6)  THEN 
ETA.O=ETA-(ETA-ETAO)/2.1 
GOTO  20 
ENDIF 

THETA=3/(TE**2)  -  2*VE/(TE**3) 
J=F*DSQRT(KT*TE/(2*PI))/(2*PD*DEXP(-PHI/KT  +  F**2* 
1  THETA/(24*KT**3)) 

J0  =  J/J1 

D=2*DSQRT(2/YE)*TE/PI 

TEST1=1/YE 

TEST2  =  1  +  F**(1./4.)*D/(PI*(D-1)) 
TEST3=-F**2/(8*KT**2*TE**2) 

TEST4  =  -PHI+KT/(1-F/(2*DSQRT(2*PHD*KT*TY)) 
IF((TEST1.GT.TEST2).AND.(TEST3.GT.TEST4))  THEN 
ECODE=0 
ELSE 

ECODE=l 

ENDIF 

RETURN 

END 
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APPENDIX  J  ARC  P80GBAM 

************************************************** 

* 

*  ARC. EXE 

* 

*  PURPOSE 

*  DETERMINES  THE  CHARACTERISTIC  VALUES  FOR 

*  THE  CATHODE  SPOT  OF  AN  ARC  BURNING  ON  A 

*  COPPER  PLATE 


DESCRIPTION 

LO 

PO 

MI 

El 

PHIO 

BETA 


OF  VARIABLES 
HEAT  OF  VAPORIZATION 
PARTIAL  PRESSURE  CONSTANT 
MASS  OF  COPPER  ATOM 
IONIZATION  POTENTIAL  ENERGY 
WORK  FUNCTION  (EV) 

FRACTION  ENERGY  REMOVED  WITH 
COULOMB  INTERACTION 
THERMAL  CONDUCTIVITY 
ERROSION  ENERGY 
AVERAGE  ION  CHARGE 
AVOGADRO'S  NUMBER 
ELECTRON  MASS 
ELECTRON  CHARGE 
BOLTZMAN'S  CONSTANT 
PLANCK'S  CONSTANT 
PERMIT IV ITY  OF  FREE  SPACE 
DEGREE  OF  IONIZATION 
VAPORIZATION  RETURN  COEFFICIENT 
ELECTRIC  FIELD  AT  CATHODE  SURFACE 
EROSION  RATE 

SCHOTTKY  CORRECTED  WORK  FUNCTION 
PI 

VAPOR  PARTIAL  PRESSURE 
RADIUS  OF  CATHODE  SPOT 
PLASMA  CONDUCTIVITY 
CATHODE  TEMPERATURE  (K) 

POTENTIAL  DROP  IN  QUASINEUTRAL 
PLASMA 

CATHODE  POTENTIAL  DROP 

TOTAL  POTENTIAL  DROP  ACROSS  ARC 

TOTAL  CURRENT  DENSITY 

ELECTRON  CURRENT  DENSITY 

ION  CURRENT  DENSITY 

CATHODE  TEMPERATURE  (JOULES) 

PLASMA  TEMPERATURE  (JOULES) 

ARC  CURRENT 


i 

! 


ft 


ooooooooooooooooo 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


ND 

DTK 

COULOG 

ECODE 


TYP 


NUMBER  OF  DATA  SETS 
CATHODE  TEMPERATURE  INCREMENT 
COULOMB  LOGARITHM 
ERROR  CODE 
0=  GOOD 
1=  ERROR 

TYPE  OF  EMISSION 
1=  THERMIONIC 
2=  INTERMEDIATE 
3=  FIELD 
4=  NONE 


SUBROUTINES  REQUIRED 

EMM(T,  F,  PHIO ,  J,  ECODE,  TYPE) 


****************  ********************************** 


REAL*8  LO ,  PO ,  MI,  El,  PHIO ,  BET  A,  XI,  EER,  Z,  NO ,  M,  E,  K,  H,  EPS 
RE AL*8  ATP, ET A, F, G,  PHI, PI,  PTK, R, S IG, TK, U, VK, VK1, VK2 , VT 
REAL*8  S1,S2,S3,S4,S5, J, JE, JI, KTK,KTP,I, DTK, COULOG 
INTEGER  ECODE. TYP, N,ND 
OPEN(6,FILE='ARC.DAT',STATUS='NEW') 

OPEN(7,FILE='PRN') 

WRITE(*,*)  'ENTER  NUMBER  OF  DATA  SETS:  ' 

READ(*,*)  ND 

WRITER,*)  'ENTER  INITIAL  PLASMA  TEMPERATURE  (EV):  ' 
READ(*,*)  KTP 

WRITE(*,*''  'ENTER  INITIAL  CATHODE  TEMPERATURE  (K):  ' 
READ(*,*)  TK 

WRITE(*,*)  'ENTER  CATHODE  TEMPERATURE  INCREMENT  (K):  ' 
READ(*,*)  DTK 

WRITE(*,*)  'ENTER  ARC  CURRENT  (A):  ' 

READ(*,*)  I 
WRITE(6,900)  ND,I 
WRITE(7,900)  ND.I 
MI=63.54 
L0=300300. 

EI=L2377D-18 
PHI0  =  4.4 
XI=40. 

BETA=3.21 

Z=L8 

N0=6.0220D23 

M=9.1095D-31 

E=L6022D-19 

K=L3807D-23 

H=6.6262D-34 

PI-3.141592654 

EPS  =  8.8542D-12 

EER=L0  /NO 

S1=-L7520D4 
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S2=-L21 
S3=0  - 

S4=13.21 

S5= 133. 322368421 
VK=15. 

MI=MI*L6726D-27 
KTP=E*KTP 
DO  20  N=1,ND 
KTK= K*TK 
JE=0. 

PTK=(10**(S1/TK  +S2*DLOG10(TK)+S3*TK+S4))*S5 
10  ATP=1/DSQRT(L  +  PTK/KTP*(2*PI*H**2/(M*KTP))**(3./2.)* 
1  DEXP(EI/KTP)) 

ETA=(8./3.  +ATP-DSQRT(5.*PI*KTK/(6.*KTP)))/ 

1  (8./3.  +ATP+DSQRT(5.*PI/6.)) 

JI=E*PTK*ETA*ATP/DSQRT(2*PI*MI*KTK) 
F=DSQRT(JI*DSQRT(8*VK*MI/E)/EPS*(1-JE/JI*DSQRT(M/MI))) 
PHI=  PHIO  *E-  DSQRT(E**3*F/(4*PI*EPS)) 

CALL  EMM(TK, F, PHIO , JE, ECODE.T YP) 

R=  DSQRT(I/(PI*(  JE+ JI))) 

VK1=((JI/JE)*(EI+2*KTP)+BETA*KTP*(L+JI/JE)-2*KTK)/E 
G=PTK*(1-ETA)/DSQRT(2*PI*MI*KTK) 
VK2=4*H*TK/(PI*R*JI)+((JE/JI)*(PHI+2*KTK)+(G/JI)*EER*E 
1  -EI-2*KTP+PHI)/E 
IF(DABS(VK2-VK1).GT.  LOD-1)  THEN 
IF(VK2.GT.VK1)  THEN 
VK=  VK1 

IF(DABS(VK2  -  VK1).GT.  100)  THEN 
KTP=KTP+100.*K 
ELSE 

KTP=KTP+DABS(VK2-VK1)*K 

ENDIF 

ELSE 

VK=  VK2 

IF(DABS(VK2- VK1).GT.  100)  THEN 
KTP=KTP-10L*K 
ELSE 

KTP=KTP- DABS(VK2  -  VK1)*K 
ENDIF 
ENDIF 

IF(VK.LT.O)  VK=DABS(VK) 

GOTO  10 
ENDIF 

VK=  ( VK1+  VK2  )/2 
J— JE+ JI 

COULOG=  23.0-  DLOG(DS  QRT(PT  K*  ATP*  ETA/K)/(KTP/K)**2) 
SIG=L53D-2*(KTP/K)**(3./2.)/COULOG 
U=(I*E*G/(2*PI*SIG*R*J))*(L+Z)/(1.+(E*G/J)*a+Z))- 
1  BETA*(KTP/E)/(1.  +(E*G/J)*(1. + Z)) 

VT= VK+U 

WRITE(6,910)  TK,KTP/E,VK,U,VT 
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WRITE(6 , 920 )  PTK AO  1D5 , ATP, ETA. JI/J. R 
WRITE(&,930)  F,  JE,  JI,  J 
WRITE(6,940)  COULOG,  SIG,  G 
WRITER, 9 10)  TK, KTP/E, VK.U.VT 
WRITE(7 , 9  20 )  PTK  A  0 1D5 ,  AT  P,  ETA.  JI/J,  R 
WRITE(7,930)  F, JE, JI, J 
WRITE(7,940)  COULOG.SIG.G 
TK=TK+DTK 
20  CONTINUE 
STOP 

900  FORMATdX I3.5X.F4.0) 

910  FORMAT(1X,F5.0,5X,F6.2,5X,F6.2,5X,F6.2,5X,F6.2) 

920  FORMAT(1XF4.1,5XF6.4,5XF6.4,5XF6.4,5X1P,D10.4) 
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